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ABSTRACT 


Lew eexisitivity ciigital are syntbesixed boiaa 

by direet and IMiraet watbeds*. SewmL ©lassaa of these 
filters «ure iPtll kfiown«; It is also obserwd that tlieir 
strm'te^es bear oloiMe sisilirity witti one another* Ware 
digital filters, lattioe digital 'filters, LBR 2^ir based 
filters axKl or^io^mal filt^rv are all oonsid«r^ and Iheir 
design prooedtires are explained with, illustrative esES^les, 

A detailed review of low sezisitivity' digital filters is 
3 ^sezitad fr€si a unified point of view in wMob the iw»ti^ 
of orthogonality plays a eentral role* the work Is aUwid at 
givisig to ttese not faidliar wl^ iiim snbdeot, an' insigiit 
into theory and design of these filters* Wo prior baokgrotind 
in analog filter the^ory is asmAed* To Maintain flow in the 
iMact, details of nathenatleal analysis are given Mily if 
essential for gaining insight* Resent work mz ViMX inplesf^nka* 
tioii of these fil'^s is also disonssed* 



Sensitlirity an iapcirtant e^nsidtration for tiie 
flosign aind |ji|>3.e»@ntatlon of digital fil.’Inrs* In Imrdiimro 
appXioationSp low sens! ti’rity to eooffieiont walnes is the 
designer*® priwary eoiasem* fills is so l>e<»us® oost of tlwi 
filter osn be rednoed by nsix^ registers of ninljiel length 
to store its parsiaeters* SlmHarly for the ¥1^1 iapleaente- 
tlonn fillMr siruotures that are modular in zkal^ire as well as 
nwaerically sl^bli! are generally preferred# It is known that 
the eonwentioml sohemMi# based im faotorisatlon and imrtial 
fraotion expansion of the transfer fisiotiim# do not yield 
stmotnrOhi with these desirable pr<^per1d.es# 

Several methods hire been stiggested in th@ |«tst for 
the design of low sensitlyitf digital filters# An Indireot 
reallaatioB sobose based on the design of doubly tenslnated 
IC filters was proposed by Fettwels [l]* This method giws 
the so ealled ware digital filters# their low sensltlwlty is 
ei^lained using the oonoept of pseudolossless [2] elements 
analsgous to the «^noept of lossless analog elesientSfl latar# 
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of ortSiogoisality ■was la'fcsr tissd Isy Itewililt st* al# to daalgii 
©asoaded orthogonal digital filters [sj* This loatliod of 
realisatioii is based on a direet faotorisa'tion of a lossless 
transfer rntrix* Ano'tber sohuie ws proposed by Taidyana’t^san 
ai»i Mitra [6] which geneimtea a oasoi^e of lossless bounded 
two-pairs eac-fcrac-lsed from the giwn transfer function* Wcflflc is ' 
also being done to inwesttgate the deslpi of low sensitivity 
sti'uo'lufe© that are suitable for 'fiSZ ii^leiMn'taiticMa of digital 
filters [7]* 

A vast amount of Htiratwr© Is available on these 
realisation scheaes. Since all of ia»se i^thods give low 
sensitivity Impleffien-tattmis of digital filters, a need is 
felt to evolve a comon theory about thea* The ideas, scat’teiTNid 
In mm literature In Mf format foras, tend to create' the' 
ii^pression ■thiit each of ■these m'l^iods gives an entirely new 
low aonsiti'vi'ty S'tani&'tfurwf* A closer look however reveals that 
they mn be te*ought together under a imif led theoretioal 
fraiaeworlc [©J* 

The ai» of this ttmais is In put '1^ various ideas togethw? 
to give an integrated approach f«nr designing low sensitivity 
digital filters* Tia entire presen'bation is of a laitorial 
nattre with profer esasples inslitied for illue^bratima* Aai effort 
is isade here ■to give ■the reader a co^tpaot yet ownplete pictwa 
of the recent trends in the theory and design ©f these filters* 
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In Section Ip me <SLXmm& the theory of doubly termiimted 
imalog IC filters and filters* which is needed to 

build the base f<»* an understanding of "mm digitaal filtersp 
described in Section 5* Sectioo 2 es|ilains the idba of lew 
i^nsitivity for aimlog fll'ters* and also includes a genewtl 
disoussion on tbe theory of.lmir sensitiyiiy’ digital filths* 

In Section 4* we e:^laiti one by one various direct realii^-** 
tion sol^mes for low sensitivity digital filters naaely t 
lattice digital filters* orthogonal digital filters and IBR 
two-ipair struotures, Later in Seetloii 5* we discuss a unified 
theory ttet explains the siailarity between all low sensitiirity 
structiares* We also discuss how the wthogonality property 
of a prototype building bloch is weful for iti stable ii^le** 
emntation by sieans of CCIIDXC technique* For salce of coe^lateness 
a biief discussion of issues related to VLSI design of .'low 
sensitivity digital filt^^ is ali^ included# This is doM in 
view of Interesting new structures are nowadays being 
investigated to take advantage of isM possibilities of parallel 
dati processing# 



SECTION OKE 


REVIEW OF AHAIDG FIL1EES 

1*1 INmOSOCflON 

In til® TOdem ei^ of «ligital processing of signals* 
iligltal filters hane replaced Itoe conirenticMaal analog filters 
in aost of tb® appUcations* Sei^^mneCL classes of digital filters 
have very Httl# in ccsmon with the analog filters* As a result* 
th® veil established theory of analog filters has been relegated 
to the background* It is felt* hcw^ver* that the theory of 
low-densitivity digital filters is difficult to mderstand 
without a prior knowledge of the class of analog filters with 
siallar prcqperties* 

Keeping this difficulty in alnd we first introduce In 
this section scns® of the basic ciaicepts related to the aimlog 
filter theory £ 93* We consider the design aspects of doubly 
terainatid lossless filtsrs that are known to, have low sensitiwiti 
with respect to the eleaent values £l€il« V© illustrate the 
ladder and lattice TOtwork realisations of' these filters £ 11 J # 
Finally* we Introduce sosoe fundamental aspects cf alcroa^ave 
flltira.i ^se filtetn have a coi^leta ^thacry of their mm-# ■ 

We consider only the class of sdcroi^ve fil-tors that can be ' 
s^thesiaed analc^ously to the dcihly terminated IC f ilte'rs £ 12 ] , 
the jpeas^ for Imluling alcroiiBve jCilters here is that they 
constitute a convenient pedag<^ic link between anal^ and 


filters* 


1*2 PEFIHIflOliS 


file most common analog filters with low serwsitivity 
properties are insertion-loss filters consisting of a 
lossless 2-port network terminated at tie output by a load 
resistance R 2 1*1)* Ite 2-port network is driven by 

voltage source ?g having a source resistance Rg* 


R,, 




+ | 



V, 



Fig, 1*1 A two-jprt aetwc«*k inserted between 
resis-live tmidJiatiaas* 


Tkm following are a mssber &t standard imrm £m 
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■ !£r&mXer Function 

transfer function H(s) of a filter is defimd as 
the ratio of an output port variable to an ii^ut port ■variable ^ 
!l3iese -variables can be either voltages or currents, * 11103 , if 
03Qly vol-tage variables are considered, then -the transfer 
function of the S-^port in Fig, (1,1) Is given by 


H(s) 




VgCs) 


traasaittiuaee and Reflectance 

A lossless 2«port network does not abswb any real 
power. Power is absorbed only by t^e source and the load 
resistances. According to the ataaiiffinsi power transfer theorem, 
^e power transoitted to the li^d has an tipper boiaid, say 
PjjiAX* transmitted power freqmmyf, *<», can 

be normalised w«»r,t» this iq[}per boimd. This defines a new 
function t(s), known as -the traasiait^aM« * It is given by 




^^iAX 


• 0 < 1 


idiere previous iK}tatl<ms are used, 

Refieetanee of a 2«port network is a function which is 
coi^Ieiaentary to the transeiittanee, . It is denoted by S*(s), 
The relati<onshit> between transmi-ttanse aiad reflectance is 



( 1 * 1 ) 


■t(») tC-») + $(s)9(-^)j • 1 

S^ju> 

Mgg function Qi Itie filter 


The filter eharacteristie is oonwuiently esqpressed by 
iti ely^ss fiyEntq.ti€>n aM another function K(s)* These 

functions are related, to t(s) by 


I,{a2) 


20 I»og(* 


|t(ja)l 


10 Log(1 ♦ |K(ja)| 2) 


The loss function is eacpressed in decibels* Further, 
if t(a) and K{s) are of the fona 


• efS} 

idaere pCs), e(s) and f(s) are polynomials in s and are 
related by the equality 


then 


e(s) eC*^) 


p(s) pC**s) + f(®) 



PCs) m fCs)/e(»), 


( 1 , 2 ) 


Itn, (1#2) i® often called the »smMm 0 

Sewiml generaliaed filter functions viz* tt* Butterwcr^ m- ■ 
the Chebyihev funetiini ©an be^ i^eaeii for KC®) so that all the 
desipi i^wcifiti^tlcais are pr«3perly mt 


The ratio of the ii3|>ut voltage to the ix^ut current 1® 
def line d, as the driving oolnt 




a 


reciprocal is ano-ttier funetioii called the driving point 

.adni tlaance f mpc tlon Yjj^(s), An aimlog networic is realised 

using 13iese imittance (impcKlance or adsiittance) fuiK^ticms )as 


it will be seen later* A necessary and sufficient condition 
for tile realisability of or is that it be positive 


real ^A||,25 # 

Ihe driving point inpedtence Z^^Cs) is also related to 
the reflectance* It is given by 


Zjjj(s) 


±? 

S ^ 


i ± e 
** f 4 e 


(1*3) 


1*3 UDim miMORK mm synthesis 

•Ehe iemittance function ob-teined fron eqn, (1*3) can 
be realised by sevezal alternative netSiods* east cimmn 
ones are tiie so called Fos^r* s and Cmvmr* » realisations (A *3) * 
The latter is also known as the ladder synthesis# To illustrate 
this saliiodi^ m now consider a sieple escaiqple of filter design* 


KCAfffLE 1^ 

Design a normlised Chebyshev filter with the following 
specifications 

^isClNand t attenuatidaa ^ 0*96 dB* Cut off frequency 

m 1 rad/sec* 

attenuatiaa ^ 30 d»* Cut off frequency 

m k jmd/sec* 


st^band t 


R2* l-OL 



Fig, 1,2 Third order Chebyshev lowpass filter of 
Example (1,1 ), 



Cut off 


Fig, 1,3 Frequency response of the lowpass filter 
of Example (1.1), 







Fig, 1,4 A symmetric lattice filter network. 



If passbatKireqtjdresients ar« exactly TOt f or 14ie 
Gliebysiiev fimction^ using eqns* (A^3) and CA*4) we obtain 
n > 2,97 (« 3) and C «• 0,5* 

Sol-wing 14ie Feld-tifeeller e«p* (1*2) -m obtain 

f m 2eP + 1*5 s 

■p » 1 

e • 2s^ + 2s^ **- 2*5 0+1 


By using eqn* (1*3) m get Zjjj(s) as 






i? + 48 + 
+ S + 1 


1 


llhen realised by the ladder wlhod as in (4*3), we get 

- 2H^ Cg » IFf » Rg • 1 0 

Ihe corresponding m-tm&ek. is shown in Fig* (1*2) and its 
frequency response is plotted in Fig* (1,3), 

1*4 SBiiiTOic LAmcE fiimm 

Lattice realisation of the filters is also wry 
emmtm^ espeoially idien rescamnt circuits are involved* A 
sysnetrio lattioe filter is rspresmted by two lossless 
inpedance functions 1 ^( 0 ) and as shown in Fig* (1*4) 

If source and the Seed resistances aie equal then 
it can be shown (A*4) that 



11 


K(s) 

anOy 

y(s) 

i^er® U and K are related Isy 

ll(s) UC--S) • 1 + KCs) KC-ii) 

Solving eqn®, (1.4) and (1*5) for and 2^^ idiich are lossless 
and hence odd functions in s, 

** ^ (1*ia) 

o o 

(1*65) 

Subscripts d* o denote even, odd parts respectively, 
fhe iopedsince functions are realised by any of previously 
nentioned methods (A *3)* We now oonsider an estaaple of latti<» 
filter synthesis* 

assign a 3*^ order Chebysbrr lattice filter with 
specifications as in exasiple 1*1 * 

Using eqn* (A*8) m obtaizi* similar to e«g« 1*1 

K(s> • ♦ „1,5, S 

tl(s) *• 2t? ♦ 2S^ tfS « ♦ 1 : 


R 


U +1 


0 o 


^a 


RC2, 


Za> 


(1.4) 


t( s) 


(2^ •» R)C2^ » R) 
R(Z^-2^> 


(1*5) 
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Separating even aiwl oid parts we get, 

r 

Kq »* 2S^ •*■ 1*5 a, Oq « 2s^ 'fr 2*5 », » 2s^ ♦ 1 

Eqns# C1*6a) and (1.6b) give 

» 2s and • 2s ♦ 2/s 

The ooiqpiet® lattice imtmtk is shewn in Fig, (1*5)* 


1-0. ' 2 H 



2 H 


Fig* 1*i tmpmm lathee tllter -©f EitaBi*!© (1,2), 
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1*5 MICROWAVja riJU'jaaKS 

For microwave applications,, we require filters -^t 
liave low sensitivity and are also suitable for hl^ frequency 
operations* In such applioations, instead of using lue^d 
elements si»li as X. and C* We use short circuited aM open 
circuited transiaission lines* In con-^st to the luo^ed 
ones, these elements are distrihuted in nature. Although 
the Immittance of these distributed elements involve espo* 
nential functions [ la], Ihey he converted into rational 
functions by tising Ihe so called Richard# transforma ti(»i * 

The design, therefore, can carried out in an equivalent 
domain* 

Richard# Transf ormatlon 

This is a transformation frcm the s-domain to the 'y- 
d(»Md% where 

^ » taiti(sT/2), idiere 

T •» 21/r m Rcanid trip delay in the transnisslim 
line of lergtii 1 
V • velocity of pn^c^tloii 

The Inverse uransforwitlon is 

,« :L iii 

■: Thu 
fieicti<»»s. 


eaqponeiitlJils are triussformed to rational 
;t I# also obssreed m#) that all the 
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of «1abillty are preserved in the ^ also^^ 

Unit Element 


In microwaw tfc«6«ry an iio^ortant eleaent is the a'-port 
element ©ailed ttis unit ele»nt^ vdiich has the following 
lamnsfer isatrix 


1 

_J 

z 

' Zm 

1 1 . 

cv 

1 

1 

1 


r 

ir 

1 

1 

? 

1 ^ 

2 ^ 

( 



It is symbolized as shorn in Fig«Cl«6)« Its 
isporlaiiioe will be disoussed later# First, we will 

©onsider s^aie design aspects of miorowawe filters« 



Fig#. 1#i A unit eleaant# 


I 


> V 




T 

V 


u f 


Short 

/ 

U £ 






2„ 





C») An Indmtw 


z ■ 2 ji. 
f 


OPf (V 
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g,omi<ieraMoiis 

The design tecliniqiies for Inoijed eleflaent filters can be 
stspai^t- forwardly carried over to tiie adcrowave filter design 
by virtue of Richard* s transformticm. In liiis method| L and 
CJ elesients of tiie former case are replaced by short circuited 
and open circtiited lines respectively in a meaner that is 
shown in Fig* 1*7(a) and 1*7(b)* IMs methodi however, has one 

drawback* In .the.. microwave ..filter d eiimu. .all ..elements can, ...be 

itctnnected only in parallel . 

To overcoTO this limitation, we can use Koroda* s identities 
to convert a series connection into a parallel one [12 ]• O^oe 
such identity is shown in Fig, (1*8)* 
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Fig* 1*8 An eirai^le of Knroda*s identity* 
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tiait EleTOBt Extraction 

One can avoid tiie me of Kwoda’s identity, idilcli is a 
coabersone process £ 13]# Ibenever a is^iint element cannot be 
ex-faracted, we extract a unit element from ttie iaoittance 
function, which is a PRF, A unit eleiaent extraction is done 
using the Richards* ja^orea (A ,6), which leaves the remainder 
again a PRF but of a lowm* degree as shown in Fig, (1,0)# 
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Fig, 1,9 Behmmtie illustratton of Richard’s theorem, 

complete procediare is best illustrated by an exwnpie, 

RVAMPTje 1 n 

Besign a microwave filter with same j^peoifloations as 
in example 1,1, Our aim is to synthesise the immittanee 
fimction in as few st<^s as possible, 

'We would . .. U^e to exixmct a capacitance .first, 
this can be done if we choose the lower signs in eqn, (1,5) to 
get 

Y(s) • 


t.,.i 

2ir ♦ s ♦ 1 




m 2s + 

^ere 

Y^(s) • 


ml 


2S^-l-S4-1 


It 


SiiKs# a ©apacitanee has already been ®:Ktra©ted, aow i#e tiaw 
to extmct a unit ele^nt* Hoting that 

Y^(1) • 3/4. 

w® extract a imit element with » 4 / 3 * leaves Y^Cs) as 

y , i8n^ t m. t , n 


*• 9/4s + YjCs) 

idler® 

YjCs) - aM Y^ci) • 3/s 

Oiice agaln<, by extractii^ a unit eleaent with • 5/3* ^ get 
Yl^(s) m 2/58 +1 

the complete filter netwexlc is shown In Fig# (1*10), 
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Fig* 1#10 Hierowav® filtwr of Bxail^le Cl#3)# 




tow SENSITIVITy OF flLimS t IWORHAOT COI^IDERAHOHS 


2.1 INmOOOGtlOIl 

Iii the previous seotlon, m familiarised ourselves 
with some of the important elasses of analog filters. Such 
filters are able to approaeimate the ideal loss characteristi© 
in an efficient manner. Besides this, we would also like to 
see that the practical circuit is such that iwall changes in 
the mlemnt paremeters lead to only small changes in the 
performaiMMt characteristi©. fwo different circuits realising 
the ideal characteristic may have different sensitivities. 

It is the designer's task t«» choose one that has least 
sensitivity* fhm sensitivity is an lapcartant consideration 
for the design of any kind of filters. 

In the discussions that follow, m first see how «je 
notion of * sensitivity is mathematically qmntified. 

This imuld hilp us understand the low sensitivity of analog 
doubly te'rmimted lossless filters in m^cise terms. We will 
also intrcKiiroe Ihe concept of Istructural* and * lossless' 
boundedness. These notions are needed to e^lain a ^ntral 

about the low sensitivity of digital filters that has 
been evolved in recent years. Finally, we touch tjpon the 
vec1»r space model of signals. latter is important today 

due to tae concept of parallel data processing idiich is now 
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possible through VLSI li^leiaentatioEB of digital filters. 

2.2 lEFIHIflQH Cr SEUSITlVIfY 

If H denotes 13ie transfer fumtion of a filter and 
p Is one of the design parameters, then the sensitivity of H 
w*r.t. p is given by 1.14 J* 



In o-toer wards. It is the ratio of fmctional change 
in value of the transfer function to that of the design 
parameter* Sensitivity is a figure of merit for any filter 
struotwre* 

In case of analog filters, the parfoxwance characteris- 
tic should not be affected mtssh by elesent valtes due to their 
tolerances* In case of digital filters, tolerance is not a 
problem because multipliers have a predetermiiied values Which 
is stored in Mmory registers* However, these registers 
can have only a finite word lengHi, This restjOlts In tte 
qmntisation of these multipliers* thus we have at ouor 
disposal only a finite nisd^er of realisable transfer func-ti'Cns 
beeaimie tlie set of paraswter values is also finite. These 
transfer functicwui shcmld have low sensitivity, if the 
multipliers are to be iiplemented with only finite precision, 

2*1 urn gMifiviCT or ymsms niasm 

It was seen earlier £l$ } that in the ease of dmd»ly 


terminatcsd lossless filters, the source delivers its maxiiau® 
availahle power (MAP) at certain frequencies, if K(s) is a 
selectiv® f motion viz, Ghebyshev, Butterwortti etc. These 
MAP frequencies, shown in Fig, (2,1) by force the low 
sensitivity in the passband region. Before we state the 
reason, we would like to Illustrate tlmt the sensitivity is 
precisely zero at these frequencies, 

fii.A 

' " , _ 


” (.Of " li.'i 

Fig, 2,1' MAP fraqismcies in the pass1:^d. 


EXAMPI£ 2,1 

Hetwojrfte of Fig, (2*2) is a 2^^ order Chebsrshev filter 
idiose response is shown in Fig, (2#1), Its MAP frequexaciea . 
are • 0 ai«i • VC2L*C)/I.y/S^ rei^ctively. Show ^sat it 
has zero sensitivity at these frtquencies, 

Transaittame t(s) is givm by 


Jim 


■|(g) 2 

e a(^)* Using eqn, CJ|#1) we can show that. 
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Fig^ 2#2 Third <jrder Chetoyshev filter of Esceusaple (2*1) 


where x is siiher 1 or C (eleieent vel»«). 

a(3«^) is giwn by 
t(») 

I Sm^ 



takiiag derivative of a iif*r*t» 1., we get 


«4&»®££^^(L%^««^^I^(31>G)+C2L -C) j 


(2 *,2) 


Siibstituitlhg for <a « aM mim eqa®, (2*1) aM (2*2) 

S^(«o) - 0 

SlBllarly substituting Iw » - » again get zero 

sensitivity, lero sensitivity with respect to C can also fee 
shsmi* 

for 

In emmmtim with tlie meS/Imr sensitlvi'^* an 
ii^im*tant conS;traint to Imi jhiqposed 011 feo'^ei digital and analog 



filters, is tJbat of boundedness# We introduce this concept 
in the following subiNBctlon before discussing the general 
viewpoint of low Bensitivity<i 


2*4 BOIlKDEDtjESS OF IRAIISFER FlINCHOM 


As an illustration, w® comMer a digital infini'te 


impulse response filter having transfer functicm as 

a 4 a * 2 auA j,. . 

hCz) - , |H(e^)| <1 


1 * h^z *♦**♦, bji® 


By structural boundedij^ss of H(e) we mean that as long as 
the TOlue. of multiplier coefficients, denoted by n^, are within 

a permitted range, H«gnitude of H(s) is bounded above by unity, 
fhis boundedness ©an also be of anoHier type. Consider for 


example, liiee following transfer function i 


HCs) 


a ♦ a-a*^ ♦ ♦ a^ a*^ 

:rrt — prr — 


Hei^, the numerator polynomial is ;jfu8t the mirror image 
of tte dencminatcr polyncnial* As a result, ^Sub magnitede of 
H(a) is alumtys unity. Therefore no matter what changes are 
made in any of tbt multiplier, it is always a bounded function. 
This is Icnown as iAm losMleiy . bo«sletoi:«| ,. This kind of 
boundedness is similar to that of lossless functions in the 
analog case, 

V# will now see how tbi bowlediiess of ttie transfer 
function can induce the low sensltiwtty property in liie filter 



2*5 LOW SEHSITIVITY OF BOtMED FlM;i?IOt3S 

Let us su|ipose that tte laagnlttide of the lareinsfer 
function HCs) is precisely uni-^ at certain freqtj^ncy for 
an inf ini 1» precision iioplenentatlon of the ©ultiplierfl* If 
now a multiplier is qiiantized* then ttm magnitude of 
transfer function can not iiKsrease* this is because tbe 
boimdedmss is * enstired* * As a Hesult, the plot of H(e) 
-versus m. is concave downwards as shown in Fig* (2*3)* 






Fig* 2*3 Illustraticm of low passband sensitiidty* 


Clearly one can see that 




'9 ffi 


(£>h 


tberefo« (2.1) « get - 0 at 


In otaier wcMCdSt boundedness of transfer fimction can 
force low sensitivity in its Ij^lenentation* In ease of 
doubly teaninated lossless filters,! immmii.t'tMme Is a b@i»i«Mi 
function* this is so because mxlms pcmer is transmitted onl- 
at More is the nunber of MAP frequencies in 

-Ihe passbandf lower will be -the sensitivity of tbe transfer 
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function* Structin*al boundeciness is referred more ctmamonly 
as s tructural passivity in analog case. 

So far we bawe considered situations where tlm»re is 
only a single input sigml (or sequence) and a corresponding 
single output signal (tnr sequence) . In digital f liter S| it 
is desirable that multi-input seqmnces too ©an be haiMiJed. 
This is in order to maxlmi^ the ‘actiwe* time of imch 
processor unit especially when we consider the VLSI Implemen* 
tation* In the following subsection w® examine this more 
general situation. 


2,6 SIGML SPACE AMD SIRUCTOAL B0Um)EDKE3S 


From now ommrdB In IMs section, we ©onsid^t* signals 
as a m^ber of n-dimensional wctcr space* Correspondingly, 
the transfer functioii is a nwpcrt tfansfer matrix* This 
n-port structure is shown in Fig, (2,^)* 



'^*ip^t stni©1»re. 


Each port (1) is jepresanted by an irgjiit sequence X|^(n) 
output sequence’ y|^Cn)* 21^ relaticaisMp ©an 




also l» looked iipon as a Umar traiasfor»ati<m T from the 
inpiitHireotor space X to the output i»c1«5r space Y, As aa 
illustration for a 2-port case^T is a 2x2 matrix as giwn 
below. 
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Alteraatiwly, t^ie irput-output relationship can also be 
characterized by a chain matrix 
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For a general irport digital filter structure, 
stnsstural boundedness is equlwlent to saying lhat ^ie 
matrix T(x) is contractile £ 16 ] i*e * 

^ - X < 0 ^ ^ C2^3) 

' ZmiSr 

iftier® superscript * t’ denotes transp«>sltlon. 

If equality sign hoMs in eqa* <2*3) for aH » then 
T(») is the twunsfer matrix of a lossless filter j , !#«♦, input 
ai»a output inner products j, are ''equal*- in such a case 

• x”*X , where t denotes transpositicei 
or, 

t |yj^(n)l ^ » E |i^Ca)l * 
i i 

•« in oth.r word*. T{Z) 1* n.gatlw .^il-dofinl 
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Thu® the lossless natore of ■&© transfer laatrlK csan he 
looked upon as the gemraHsed notion of lossless boundedness* 
We will see in a later section hm tMs can proiride a low 
sensitivity and stable design for digital filters for 
parallel processing of signals* 



mcTion mm 


mm DIGITAL FILTERS t imiMECT REAIISAflOH (M WV 

sMmiTinTr filiers 

3*1 iHmoDucnoN 

We are by now aqi^inted with a olass of analog filters 
viz* dowbly terminated lossless filters* It was shown in the 
previous section that a lossless netwi»’k inserted betieeen two 
resistive terminations possessed low sensitivity with respect 
to its element^values* Ihe prcperty was seen to be -toe result 
of the striKJtural passivity of the filter neturorlc and the 
restrictions placed iq?oa it by the laaximum power transfer 
theorem* By analogy, it was later seen that the low- 
sensitivity also extended to digital filters Whose filter 
fmictions were structta?ally boin^edt 

In this section^ m discuss an ij^ertant class of 
digital filters that imitate analog doubly terminated lossless 
filters in the digital d<miin* Itie:se filters are kmm as 
•Wav® digital filters* (WBPs) and were first introduced by 
A* Fettwels [ 1 ]• Owlr® to their pcpularityi now they have 
also been included in some tesrb^iwsoks £ 17 ]« A brief overview 
and details of the current state of art on WIFs can be 
recent '.featorial paper 





In tiMS discussions, we cower ■Kfee Msic concepts and 
design techniques of these filters# Their low sensitiwlty 
pr<i53€rty has also been eatplained# Some simulation results 
hawe been included to glwe a more concrete picture of WiW& 

1o -the reader^ A^ppe^div iS> 

3#2 tmiMimmm 

Our purpose is to design digital filters that eahibit 
wery low sensitivity to clsminges in multipliers. An attractive 
apprcach is to consider low sensitivity analog filters and 
to try to find analogous structuz^s in tiie digital domain# 

If a suitable transformation can be foinwl, then we can 
realise low sensitivity digital filters by an indireot 
methCKl# This would also obviate the need to develop entirely 
2 ^ design teclsiiques in digital domain* WIFsg, however are 
still very attractive l^cause of their simple design techniqcms* 

To synthesize wave digital filters.^ a transfonwiM<m 
is required by which sue analog nitworU is converted into an 
equivalent digital structmr#' tb& treouifc»:*nation should be 
such that mmm of the stability conditioiis are violated# The 
design of ICEPs is carried cut In the transformed dcmainf 
deim^ted by 'Vp * 1%i8 domain is ^ealogbiis to iSm actual od^len 
fpeouenev Sn«dcmainA whicli is hencefcrtli mailed the referezice 
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doEMiin* The relationship between ip t # and z is gi-ven by 

'f m^lh 

or, 

M -1 

. fanh (sT/2) C3.1b) 

1 + z 

Eqn* (3*1b) is siBilLar to one used in the design of 
ssiorowate filters to conwrt expmientials Into i^itionals. Here 
T denotes a imit of delay, most often the Nyqnist saag^ling 
period. Fig, (3*1) shows that the area inside unit circle In 
tte z-^lane is aapped onto LI# in the op -plane * Therefore, 
stability conditions in the i-df»ain are preiMfrwed since all 
poles inside tlie unit circle z-Tplan® are traneforsied to LHP 
poles in -the reference dooain, a 



J^ig# 3*1 Mapping of »-plane csa plane. 



In ^ reference doaain, all porta are characterized 
by relatioiiships between ctnrrei^tonding woltage (V) and cDrreiit Cl) 
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slgmls. In the case of wave digital filters^ these signals 
are r€^resen1»d by the so called ’wave quantities* * fhese 
are the incident and reflected wave signalsi^ denoted by A and 
B respectively. Also, each port is associated with a port 
resistance R which can be arbitrarily chosen* Relationships 
between these signals are as follows 

A » V ♦ RI, B » ? - R1 (5*2a,b) 

Represen.ta.tion...of . elenenti 

t 

tim desigGi of bUFs, eleaants vis* resistors, 
capacitors, indnctors and voltage sources are treated as 
single ports* Unit elements are, however 2~portB* Using 
eqns, (5,2), all these can be trensfcnwd into 

equivalent digiial sub-structures as shown in Fig, (5,2)* For 
illUEstratiim of this, an iMmtmr is considered here, 

mmmjjd 

Find the i^uivalent ifXF repreis^tatifm of an inducte*. 

Fear an inductor, vol'toge-cmrent relationship in ^ -dcwain 
is given by ■ 

¥ m- ( I. being tte value of inductance. 

If we choose port resistance E «• S*, '^len eqja* C3#2) gives 

A m (1 11/) RI S *■ *41 R1 
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Emm B • - iX-zll A 

(1 4 */^) 

and §qn« (3*1a) giva® S ** -«"***A# ■ 

Thus an inductor in th® ^-d<»»in is aquiwlently raprasented 
hy a delay in the z«ii<min as shown in Fig, (3#f)« 

, + V - 

■WST)! ^ — A -V 0—1 1 j — B 

Fig, 3*3 IfKiuetor and its IDF equivalent of Ixan^le C3*1) 

¥e have seen hoe it is possible to simulate in the 
digital d<»ain the m&in elei^nt® of the reference filter. 

To complete ^le equiimlenee, %re also need to satisfy varioias 
•topological rules such as the Kirchoff*s lavs in the reference 
d(»ain. Similarly, •various ports in nave digital filters can 
be connec-ted either in series on in parallel. However, imless 
pcu’t resistances are identical, ttmm ports cannot be connected 
directly vithcmt violating eqn, (5,2), To overcoaws -tMe 
difficulty, a multiport eloawnt ©ailed an * adaptor* is used 
for such Interccasiections* Adaptci^ play an insportant role in 
tbe realisa-tioii of li^fs as it viH be seen in the folliwiiii^ 
discuilsion,' 

M.stat^ earlier, adaptors ma used to ccnuwct ^digital 
ports eito»r In series or in parallel £l#], 


are 



3 ^ 

basically ’standardized.* digitol structyres that consist of 
only mnltipliers aM adders. There can be both serial and 
parallel adaptors depending upon their ftinctions* These 
standardized n-port adaptors are shown in Fig. (3*#) * To 
illustrate how a typical adaptor equations are obtained^ we 
consider the following exau^le* 


Design a parallel adaptca* for interconnecting two 
ports with port resistances and R 2 res^ctiwely, 

Eqns. (3*2) glm& the general deseription as 




®1 • ^1 

- I,R., 


^2 * ^2 **" ^2^2* 

®2 • ”^2 

-1^2 

For parallel interconnection of pcrtSt we imim 


* ^2 

1^ w — 1^ 



m * B2 



amdg 

Am ** B* ** B^ 

*^2 

m D ' 


Solving eqn, (3*3) foP and w® get 




% 



(3*3a) 

{ 3 »»>) 


Ra 



(a) Unconstrained parallel and series adaptors 



(b) Constrained parallel and series adaptors 


3*4 Standardized N—port adaptors. 
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2G- 

^1* - mpny 


*» l/R^f Gg • l/Rg 


the 2-port adaptor is shorn in Fig, (3«§) • 

Ai) -(£) i 4 - A? 




V 




1 :. 


-G-~ 


Fig, 3*5 Two-ipart parallai adaptor of Example (3,2), 

Serial and parallel adapters* maed Fi aad S2 
respectiwly in Pig* (3^-) ai^ taaco^traimd > TM name is 
giiren heeause i^re is h®' restrietioB oa.the valoes of the . 
port resistances. These adaptors require ir*1 anltipllers* ■' 
each h.^ (i • 1 to n-1) feeing an indi#in^ent quantity. 

Equations in Pig, (3#%fe)» hoeeeer* suggest that i^ i» 
unity for a port i* then its corresponding refleciad sigiml 
is ind^epeiaient of the inciiimt sipial Thus the port 
feecomes VareflecMon free * , Am adap'ter hairlBg a reflection free 
port is lulled a constrained adapter and is deno'iM fey Ft m ii .* 
A censtrained adaptor requires oniy n*€ i^tipliirs fmr it» 
r«alisati^* / ' . ■ 

A reflection free port oan aleays fee coinmcted to 
any ottinr pin*t without pesstlting in a ''delay free loop. Such 
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l©«^a should be absent to a oc»^utable digital struelaire. 
Since there is no path frc®i the iaput signal to the output 
signal of a reflection free port,, a delay free loops in this 
case are autoiBatically airoided« If certain precautions 
(A ^7) are kept in Bilndg then ci^utable Wf striK^tures can 
be obtained. 


Baling considered imrious building blocks of a WW 

Structure, we now consider how a tcm^sfer function denoted 

, / ¥ 

by Hj|( ijn ) in the reference dc®»ih/rwilised in the digital 

dmain, Sisilar to the analog casein r^jialisation is 

Imsed on the ^sign of a 'psutdclossless* network teraina1»d 

by an ‘ equiwaleat* resistlwe element (<^n circuit in the 

im design). 

Consider again the netisnk: of Fig, the ir^ut^ 

output port equations mm be written as 

*2 • »2 ♦ • Bj . Tj - 

iinct'liie output port is resistiifely teradnated, we hare - 

Wtst'^maTf 

Aj| « 






( 3 #^) 


Xt Biz) denotes the transfer function in digital 4<»ln then 
eqn* <3*4) gives 


Biz) 


B- ? 

. ^2 I 


'ArO 



m 2Ha( 'f )! 

2 - 1 )/ ( 24 - 1 ) 


Ihe 2-s>ort digital structure is In Pig, (3<ii) which 

realizes the transfer function H(z), 




(S, — <- 


Di6*i tal 
nvo 

Port 








Pig* 3*6 A two-ftort digital struntw^* 

3#S L0W BEHSmyiTY QF 

!lhat wave digital filters show low ses^lMvitsr to thi 
e^iffici«ats* i|uantlj^tion can he intuitiwly understood, 
fhis is so because IDfs are, oblained by a bilinear trans^ 
fcnation of the reference filter leaving this property, 

to ea^lalja the low sensitivity of llBPs iian:*e rlgoroiwlyf 
the concepts of*pseudolosslesin^ss* and * pseudopower* need be 
defined for wave tuantities in ^digi'til doaaln [ a ].« 

Consiiter for a 2-port|f the total nsuedopowar P, 
absinrbed by it iddch is- defined as 

F • ■ C 3 #S) 

where E is the port ihisistanae and a and b denote input icii 
..ou%ut wave qutmtities respectiv«ly« In the case of 



doubly tenaiimted filters^ it was soon In Section II that 
the inserted network consisted of lossless elements alone, 
fliis resulted in HAP frequencies in j^ssband wltb zero 
sensitivity at these points* On similar growids it can be 

shown by an identical reasoning that if the elaaents in WF 
stmctoes are pseudolossless tJ^n «DFs have the low sensitivi-ty 
property* ¥e consider as an illustration, ihe case of a 
i**port series adaptor* 

mmm . ,?«,? 


Show that a 2-^ort series adaptor, as shown in Fig* (3*f), 
is pseudolossless, lue* it does not absorb any pseudopcn^r. 

A. ^ p- -<4) ^ A,, • 

'f . y ' 

I • 

B, — — i 0 

fig* 3*? 'fwo«port 0^*108 adapts of Exmaple (3*3)* 
the pc®’t equatioms are given by 


• (1 ^ k) 

^2 » ^2^) (1 •A)Aj^ (3*S) 


fetal power absorbed by the 


is 


S 2 

ii,,j 

s 








(f*T) 


in iSkf) and siaiplifying* we get 
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P » C(2*^)A^ — XAg) (|^ «• 

Siibstitiitlng for A in the last tera^ m find that P • 0, 
fhw the adaptin* is a pstiedoiossless eleaeat. 

Similarly# other elements incltiding delays too oan 
he shown to he pseudolosaless, A rmiator^ i#e# an open 
oircuit however# is seen to be a passive element. Therefore# 
idU. results of the analog ease ' apply here' as well# ' Henee 
Mse low sensitivity nature of wave digital filters is seen 
to he preserved. 



W# n^ eoiviider various kinds of VW structures Idsat 
are based 4^ . the synthesis of analogous reference dUseain 
filters, ^ These filters were disoimsed^in Becti^cai I ^and ware 
later shown to possess low sensitivity property, liaaely# 
we can design ladder or lattice wave digital filters. Besides 
Ittase# unit elment based lliP are derived froa the 

ocrrespcaiding microwave filtai^. 



Ladders are the most tmmm stnsstures imt wave 


digilal fil'^s. In these s^g^tures# serial ant perallal 
adaptors alternate# imitatiag tiie series and shunt eonneotlons 
the reference filt^ as shm^'ln Pig* laolt adapter 
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is a 3-i>ort, the second pcsrt representing either a 
capacitor or an Inductor* In MDF^ an L or C element 
isvolT»s a delay, Thm tim structure is canonic in delay 
eleEsents* The 3rd port resislsinc® is so chosen (A ,8) that 
each adaptor except the last one is a constrained adap'tor*£2©]* 
This imkes the design of adaptors possihle hy using only one 
multiplier each* Since tte last adaptcsr is unconstralntdi 
total number of multipliers needed are nel* To lllt®trate, 
u® consider an exsnple of ladder WW s'hructure. 




Fig.* 3*8 A ladder based yUM stna$ture, 


/ The filter netwmrlc cC Fig* (1*2) can be redrami as 
Fig* (3*9)* Derive Its equivalent USP ladder structure. 



^Ig# 3*9 Third order Chsbyalieir filter of &eaaiple (3*4). • 





flue equiimleiit WDF striKjtur® ia aaaily drawn as in 
Fig, <3,10), 
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Fig* 3*10 Ladder l:«ised WF »trmt\we ©f Eswa^le (3*4) • 

C!3©osiiig iK>rt 2 as 'th® deijendeat p«rt for sash adaptor* lie hS'W 
Adaptcn* 1 I i^pe SI 

^ § ^2 * ^ * ^3 * ^ 

^11 * 

Adaptor 2 i type Pi 

CSi.| • 1/3 § Oig •» 1 'n Cl*^_ *• ■ 4/3 

>■21 “ 

AdaptPf 3 I type S2 

■R^ *• ^4 ^ « 1 ■ 'g, R^ •• 2 

^31 * iz/s * ^3® * 


It is observed that tlie tliird adaptor requires 2 
multipliers because it is uDco!]»traiued as none of its port 
resistances can be arbilrarily chosen, 

toit ele^nt based structures 

An attractive structore fcxr mve digital filters 
is obtained by the correspontiing microimve filter design 
[ 21 3 * Since unit elements In mm digital filters are 
i^alised by half delays^ these delay based structures have 
several advantages, Firstiyjp the delay^free loops which may 
be gemrated in otbar desl^ss are very easily avoided. 
Secondly^ use of delay elements reduces the number of 
multipliers in each delay^fTee path. This saves the ocM^uta^ 
ti«m time reqtilmd fee* the promsslng of each iiqEJUt data. 

In the followiz^ ejnple we comider design of a unit 
elenent based 

mmm . . 

Design e unit elmmt basi^ iDf trm the reference 
filter ®f Fig, (1,10), 

Referring to Fig, (1,10), its equivalent Wi struclw® 
is drawn in Fig, (3,11), 

In this case, it is observed that all adaptors are 
woonstrained (type F2) , Opecificaldoias of these are ' 







Fig* 5*11 ltei1?-eleaient baaed »tnieto« of Ejcaiapla (5*5)* 
Adapt^a* 1 * 

- I, G^ «• 2, G^ « 3/4 
• @/i5, • a/5 

Adaptor 2 I 

•» 3/4* Gg • f/4^ G3 i. 3/5 
«. 5/12, » 1/3 

Adaptor 3 I 

'G^ m 3/5, Gg • 2/3, G^ «* 1 

' • 9/17* X^ •. 1©/17i 

A point to nolo, aa sboiiii in Fig* (3*12) ia tint 
m&f two lialf dala^rs preoeding a digital . port oan bo ddiood 
'tonotboi* to f om a fall dola? olMMitit* TbxM in aotyal 

; - .. ;/;.' ’ . ' ■ ' ■ ; ' : : 




design, no half delays are lased^ digital struotwre of 
j^ig# (3*11) caa# tlieirefCTO, be as Fig, (3,13)« 



1 1 

I 

“ 

N 1 5 i 

N 1 

---4 

! 1 

1 <-[fi<-l 
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f; 
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Fig» 3 #12 Equivalene® of half-delays aixi a full-delay 
element. 








-> Bs 
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J^ig# 3*13 IJalay based ii©F stx*u 0 ture of Ixaiple (3*5) • 


DlEital, Filters 

lattioe waim digital filters are i^is® aa isp^rteit 
olass &t eew digital filter#* Ohtaiued from the Mttlee 
rafereno# filtxwi in similar famMott as al^ara* ^y heooae 
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fmmr multipliers in ttie digital domain t22 I* 

Consider lattice network of Fig, (3,14), 



fig, 5,14 A syBTOtrio lattice network. 

Here both and are lossless reactance functions. 

Just as we had used earlier an iapedaiwe transfer matrix 1,23 
t&t Idle mmlm case, her® we define a scattering laatrix [s] 
which relates ii^ut and output wawe Qiantlties by 

I - i 

where, 

I • t®1 ^ 


and 

‘ r* 

®11 ®12 

[S] • 

®21 *22 

It is observed (4,9) that for a sywisetric case 





Hence I the transfer function is giwn hj 
H(z) - — a^. - « - I 


^ A *0 " 


The timnsfer function is realised as in Fig, (3*15) 


A 1 ^ 


0()i. 


h 


Sb 




5‘ig* 3*15 A lattice realisation of wave digital filter 
Both Sg^ and s^j are ohialned hy the transf onaation 
of reactance functions* Consider* an exaaple for illustration* 


EXAMPljB 3.6 

Design a lattice wave digital filter for the specifica- 
tions as in exai^le (1*2)* 

In this cascif and are obtained as 

m 2s Z^ It 2s # 2/s 

After transforraatton* they are redjmwn in Fig, (3, IS), 



^ 5b ' ^ ^ — J ' ■ 

Fig* 3*16 Lattice IfflF structure of Exasqple (3*6), 
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Adaptors in Pig* C3*,16) haw; tlse following par*aBieters 
Adaptor 1 I typi« P1 

« 1 , • 1/2 and A-J 2 » 2/5 

Adaptiar 2 * type S2 

^<1 •1 ji' • 1/2 y Rj 2 

^21 * ^>22 » 2 / t * 

A« a wariation, one can also apply Ricliard* s tliaoi^a for 
Realizing and which are loaslaaa ftaictions* nils gives 
a general a^thod of designing lattiee WEPs by using only the 
unit ela^nts*. < ' 


a.ta.te. 3mm.... mm.... m 


rnSmSSSmmk 


Wf stru^tixres disomsed so far ymv% tesed on direct 
iiipleiainlation of transformed structures in the digital 
domain* State space were digital fil1nrS|, howev^^ are 
ia^lemented using Idle sts'^ variables £25 ]« 

In general state space digital filters are expressed 
by. Idle following difference «tiiatiOT»# ■ 

£(i»l) [f J Kn) ♦ [Bl 2<“) 

•• [C] £(b) * Id] i^b) (3.9) 


liliere J aiMi U are the slmte* emtput and input vect<a*s 
respeetiwly. 


A DC®’ stnjKsture derived ireaa the referenae filters 
by any of the earlier ae-^Kids oan be loelced t^on as an (11+2) *• 


digital multiport. Here W Is the total number of delay 
elements correspondix^ to i:Sm Indisetisre and cai^oi'tors* If 
these delays along with the input and oui^ut p^nrts are 
taken out as shown in Fig. C5*1f) thKi Ihe entir# network 
©an be considered a g&mmmlXmd (H+2)'*port adaptxn*. 


^ ‘-L. ^ 


1 ^ — Y 


A I 


n 


-| C M C R A L I ■G E T) 


fKt-f2)' P£)RT a "DAP tor 


!• — ^ B : 


A 2 


4 , — y 
H 7 r"" 


"IP 


fig, 3#17 Generalised e^pixi3P for 


Fen* a port i^ tersdnati^ In a delay el«men% the c^nrrtipondii^ 
equatlxm is 


or. 


A|^ '•» z 


Ai^Cw^l) • ^») 



If we ohoose A^,* s (i • 1 to H) as H state variables, A^ as 
ii^ut and Bg -as output* then the witire W- strueture ©an be 



as a 




this is iXliaitrated by 






Rualise the ¥DF ©f Pig. (1,2) as a stat» space digital 
filter (SSDF), 

WF to be modelled as a is weHlraim la Fig* (3*18) 


A:, 


} 7 


, l! A.-i 
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1 

1 , <■■■ • 

1 

{ ’ 

0 

* A.;. ; 


4- 


A ,, 


A) 




( 2 ) 




A - *■, " /'^ *■: 


Fig* 3*18 WF stmotunl for SSDF realisatloa ©f Eatai^l© (3«7) 


If ij^ut and outpiit nmVB quantities ©re denoted by 
and respeetively* where 

i «• adapter zuifibi^ 

5 • |iwt laiabw?* 

Then fro® Fig* (3*18) it is seen that^ 

®13 * Si* S 3 * Si* ^23 ®3i* ®23 “ Si 

Replaeing ^33 ^ ‘^33 elimiiiat e 

the mg&tlw sign befwe dela^,,we get 
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/Vb. O 3 2 0 


fil»™ 

33 

■ -|*22* 

A *11 - 4 *i2 

4. Jr A 

M 33 


®32 

■ A*33~‘ 

i *^22 * 

^*11 ~ 

4*12 

(3.11) 


Using tl for [X^ Xg X^]”*^ for mm Y fat 

we wrLtot the state equations frm oqm* (3#9) and (3.11) as 


|)c:i(»i)|; 

1/15 

-4/5 

-4/15 ^ 

X.|(n) 

' ) 

i ■; -46/151 

.Ej<DH)i 

' ' 1 

- 2/5 

1/5 

■ 2/5 1 

! X,(n)j 

i ' 1 

1* -t-l-^/s 1 

i: 

X-jCml) 1 

-4/15 

•4/8 

1/15: 

' 1 

1 XjCn)! 

1 1 1 

1 1 V15 1 


*1«0 “ t-^ -f As jx^Cn);! ♦ t^S (5.1a) 

IXjCo)! 


Wlieh is 1^ desired state spase fcraiiiation. 

In this seetloiig m fectissed avot attention on an 
liBportsimt Glass of low sensitivity digital filters viz* wave 
digital lil^s* It m» shown how a large , ntsaber of possihle ■ 
dealgiis in the referent doMtim are also nsefnl for the^ design . 
of digital filter stmastares' using wave quantities* *I!h# aethod^ : 
alttatgtm indireoti also gives good iwrfomaiioe with respeet 
to (i^nantizatioii of wiltlpliirs» In next seotitm* we will 
ooMldar sows ^direot approaohes to Ihs design of low sensitivitf 
digital tutors that do '' not. require a prior kiwiwledge of analog 



SEGtlOK FOR 


PlBEgf REAUSATIOII OF WV SEIgm^ITY DIGim nUTERS 
4,1 imODtICTIOlI 

In the in^idLons e#etia% m diB&mB&d the thec^y 
and design te<^inlqiies of m'm digital flltefs* the synthesis 
of ^bese filters i^nlres a kncmledge ©f analog filter timmy 
as Mil, The nethod of reallsatlfm is an indlreot one heoanse 
any s-domain transfer function^ HCa)||f is first transfcaciBed 
Into its oorrespcmdlng transfer fnnotloii In the reference 
dooaln. It Is then realised ^ ..in reference doeain and 
transformed back into the 38-*d©«aln# It nay sc»etijB^s happen 
that no low sensitivity referent fll^ exists for a 
specified hC*), Tlais the indirect i^thod is jaot universally 
applicsahle. Besides thiS|^ the <^ioe of digital filter 
struotures is liaited by the nmaber of different reference 
filter networks for a given traiafer fuacticin. Further^ in 
this apicoaoh the optiwizatjcn of nullipliersy adders and 
delays etc, is oarried out only for the building blocks vIe* 
adaptors etc*, and i«>t for the mtire filter stnasture as a 
whole^ The sethod, therefore^ gives only a suboptinal reall- 
sation for low sensitivity digital filter ©haraeteristios, . 

In s«©tioBj|. we ermine sene of the direct 
realisation awitbods for low sensitivity digital filters. We 
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start with a ladder/latticc stnictara proposed by Gray aind 
Markel £ 3 ]# We coiisider the concept of 'eabedding* of 
a boiwled function into a siiattering matrix for its 2-port 
structure realisation. We also discuss how this natrix can be 
syntl»si 2 ed by the so called ISl-extraction approach £ ^ j * 
Later in tdie discussions ^ we coKsider the realisation of these 
filters by cascading of *or1hftgoi«l* blochs timt constitute 
what are known as orthogonal digital filters £ 5 J* Finally 
we investigate the coMWtn featares of all of these TOthods and 
then try to establish a liiA: between them by presenting a 
imified approach for the design of low sensitivil^ digital 
filters in the next section. 


4 f% t fWt tAWfitI 


Pirect realisatiini of low senBitiifi% digital filters 
was first attempted by Gray and Markel £ 3 ]* An algorithoi 
proposed by thee giiras a c^soaded dig! til. filter structureii 
Each building block ©an be n^lteed either as a * ladder* ‘ or 
as a *latti@e* structure in the digital dwain* to illustrate 


this, consider a transfer functicm H(a) giiren by 

% 

hCs) • 


(4*1) 


where, 


t * 


W 

mf 


M' 

m. 


m 

t 


p. 




m e 1-2; 






v*> 


M 

1 4 2 

nw1 


ft 



transfer fimettoia is iljypaetlf raalisalQ.# by losing 2M+1: 
multiplier coefficients but tbis <toes not yield a low sensi- 
tivity structuoro* As we sbtll presently seei^ GRAY 
mmsVs algorithm also requires »l4l parameters aM is, in 
that sense, canonical* let us defim intensediate polynomials 
Aj^Cs) and B^(z) of i3m forms given in (4*1) for s » 0^1,2,** *15* 
Then the algorithm (A, 10) gives the following realisation of 
H(z)* 

n zBiz) . 

“ A 'Vicfn ^ 


Here are called 1*® tan paramotBrs # 
stage m are related (An-ll) as 


^®4l 


^ \ 




m 




®«»4t 


— 1» **1 

s kjg s 

■j 

1 

» 


1^ are called teBSSSaiMi# ^ algorithm (A *10) gives the 
values of tap parametiPi and k-parwters for any bowided 
treiMifar fmotiimk* If XCs) awi Y(s) denote. ii^ut ai^ 
ou'^ut sequencee respectively then YCz) can fee written as 

Y(*) • aCs) xC«> ^ ^ 


C4#5) can be realise# d» « 
is shown in tig* 


.ce striictiar^ irtilcli 
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Fig# ^*1 fiKi-lwt iattie® sianiciair®# 

Using tMs lattice structure ani. eqn* (4,2), «e can 
realise tbe cascaded fora of digital filter structure as siiown 

In Fig, (4,2) , 

^ i 


‘M-l 


65 m ; 


<• 


(5i 


jii 


H 










c}-- 


'j- - 

/k, i 


1 _ 

G'to 


r ' 

\ 












^ Y 


■ Fig, 4,2 Caseadfid digital filtinr strucfer® of Gray and Marfcel, 
Aitm*natl’¥®ly, a rearrangefflant of teros leads to a 

ladder »tn»ctwM (fIgDM (4,3)). ^ transfer aatrlx of a 

2-port ladder is giwn as 



I 

a 1 

<*> I 

*®at‘1 i 

! 


'n 


i \ 




2 


(4,5) 


I B, 


a 


A...y 

A(\t 





---0 


kiV 

A M 




\J 

vi 


A-i 


2 l3n,,..X ^ 


Am 


-<Y— 




A M 


4,3 f'iifO*''l?M3rt l a d d MT struct®*®. 
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By u«lng the concepts of * pseudopower* , sinilar to that 
in case of the waire digital filters, it can be shown Idbet the 
GRAY and KARKEL stnssture is also stnie^irally passive and 
therefore, stw>ws low sensitivity t© changes in coefficients 
due to quantisation [2^ 3 * It will be Illustrated later in 
this section that this is due to the • orthogonality* of the 
transfer Mitrix of the ladder or lattice building blocks* 

Let us now ccmsider an eacax^le of this al^orithn* 


EXAMBU; 4 1 

Design a Gray eund Karkel lattice filter witaa following 
transfer function 

¥imm^ 

pg • ,4 ♦ ,0 ^ 

Frctt steps (1,2^ 1#3 and 1^4) algorithmgives 
*®2 • 

lt| ** ' : 

^at “ %2 

Steps Cl»») safi / 


H(s) 


„Mii 

(1 + 0#89 n * ♦ 0,16 z * 
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* (^2 ** « k |) 

* 1 + 0*7672415 z"^ 

* ^2 ** ^2 

■• 0*336 ♦ 0^444 

zB^ » 0,7672413 + ^ 
x*^^atizig list pns^ess, m obtaiji 
\ • 0,7672413 
» 0*444 


the lattle# stmetiffir is |jip2.®etzite<l as Jjq Pig, (4*4) « 



Aa alteraattf® awiaicid ta.rmlim a tranaXar ftmotlmi 
i« t® ♦eiBbad* it ia a digital. a*-t©3t*t traaaXear watrlai* It ; 
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wa» illustrated in section II how a lossless hoiKKied fiincti<m 
is insensitive to parang tesr perturhations# In this ease, 
low sensitivity is *enfi»’ced% On siailar grounds^ we 
oan in'teitively argue that a lossless jaatrix t 00 | which 
is a gemralized notion of lossless bounded functions, shows 
low sensitivity to its parai»ters# 


Evidently, to achieve this. It should be possible to 
embed a lossy transfer fimctlcm^ i,e, om which is not 
lossless bounded, into a lossless smtrlx# If H< 2 ) is any 
causal and strictly stable rational transfer function which 
is bounded above by unity, then it oan always be embedded 
into a matrix TCz) nMch is lossless. This result is based 
on the classical embedding procedure given by Belevitch for 
srattering fimcti<ms into a para unitary matrix [ 25 ] , 

If the transfer functlcm HCz) is given by 
hCz) • f(s)/g(E), and is embedded as (figure 4,5a) in 
the matrix TCt), then the T(z) Is given by 


idiire 


T(x) • 


fnCx) I 



fCs) hC») 

1%C«) I • I*® hCS*"*) 


N is degree of the polynomial gCs) , 


(4,6) 


holii 


For T(e) to be lossless, toUmdm relationship 
between f^g-and h 

ff, 4* lih|i •» ggt 
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Fim the knovlefige of f eiwi one can obtain fe 
using «<jn# (4,^7)# Sine# hh|^ has reciproeal roots, only 
those roots are included In h(s) that nake It a hurvits 
polynomial. For the z^^in, iihia implies that all the 
roots ar® strictly inside unit iircle, we now consider 
an example of embedding, 

4.2 

The followit% transfer functicei MCs) is to be embedded 
in a lossless matrix 

f](z) H Sia&ili. I II 

(1 4 0,8 0,25 

g t mrm giwen as 

g • 1 ♦ 0,® s*® gn • 0,25 0,8 

i m 0,5(1 4* 1,28 s"^4S*^) f* • f (due to syimmetry) 
(4,7) gives 

» gg» ** f 

• 4- i"*®) 

» t,75 B^ (#^ ♦ ®*^)]t»7S Cl ♦ ,64 B^h} 
mmB ik la given liy 

h e ,75*'^ (0*64 e 
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Hie twoMfei* mtrixp t^ierefcrep le obl^ijied m 

*'*7SCl*»‘*64 z^) *5(1+1 *28se"*^4**® 

#5C 1+1 ,282***'+^*® »75z*** (0,64+z*^ 

It can be wrified that 7^z) 7(«*’*) m 7ht» t(z) ie lossless 
matrix* 

As a 'variation HCk) can also be embedded as as 
slwwn in Fig, (4*5b) , 


T(z) 


( 1+.8 


"tsf 


+ , 25 ?^) 



(®) Esbeddl^ (b) Embedding 

yig* 4,0 Embedding of transfer fnnctieii* 

After a loss^f transfer fmxstion is embedded* into a 
lossless Mtrix* we shoidd find means of in^lementlng it as 
a digital 2«iiort udaose owizall Ascription ' is as sboim in 
Figs, (4*5)* this is what we next discuss, 

9A8CABB BIFIEimifgATOK I^..Xf^lES3.l!ATO 

Th$ lossless matrix obtained after embedding has the 
order which la same as ^ Agree of the embedAd -^Afer 
function, Qm way, thsreforei^ is to raaliA ISiis fCs) by 
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This metiKKi alas at realising anty stmeturally 
or lossless liioijnded transfer frnietion as a easoade 
of IM. 2-|)air8* ¥e start idlto a tramfer fimotion U^iz) and 
extract from it a m 2-^lr in such a mnner that the 
remaining transfer fimctltm fa of a lower degree. 

This is shown in Fig, (%#7) , 



X , 


Fig. 4,7 ExtpaeHon of stn ISR two-pair* 


If to 2«pair is ieserihed hy its chain parawaters^ 
ton to following relatlonshlf holds hetween J^(z) and 



(4,@a|,h) 

Hext, we deserihe to tmm of first and second order 
ISk 2-pair8 which are disciaan^ in detail in £ $ ] . 



The siaplest for»' hf 2*pair which can be stM^sessiwl: 
^neietted frcai giwen tspansfeir' fmoottcsi Is of first order# Ibe . 
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Chain and transfer parameters of this 2-*pair are 

A m 

i * 0 z"*^ 

HI 

mfA 

ft 

1 

m 

B «* 

+ (1— o)z**'* 

*<12 ** *2l 

C » 

+ (1 • o) 

*22 ■* * (1 ’• o)z"’'^ 

D m 

'‘•e'l 

er ♦ z 


Common dencMinator 

Cowon denominator 

mj a 

(1 + «"^) 

•• 1 + 0 sT^ (4,^ib) 


A |K»iiit to no^ here is t!iat tJit 2**pair dosoribed by 
pmrmmtmr o, has a traii«aissi«ai aero^ i,e* f ^ » 0^ at 
a • -1* Its isgjK»rla»o« idll b© showt later after we also 
describe a seoood carder SH^aJr wMeh has transiissiozx zeroes 
at z*^ « ♦ 

Order, im .2-<f*air. 

When trsnsBiissicaa mmm in ©oa^ugate pairs say at 
• • then a 2-|)air @f seeond ord^ is extraoted* It 

has twsi paraaetcrs o ai«l § with ^ Cos The ©hain ant 

transfer aatrlx has following paraneters 

A ■» i e ^Ci ^ ojz***^ e 0 %^ ** ^ 

1 *» J;C1 0)**^C|t«'***) ' % • • /0(l+2#Z*\z"®) 

C » ^ «)C‘t ♦ ■ ' *22 • "* o)z****(#^z^’^) 

!3 * o ♦ s(1 ♦ ol**^ e 



6^1 

CcMion d^nomina'fcQK* Ccraroon denwlnatc^ 

*• >/ ® ♦ 20Z ^ 4- nT^) i* t 4. 0(‘i‘t'O)z*^’KJZ^ 

(4*10agb) 

*^*B^aaval irm tiie fmmt&p FtmctS^e^ 


To imdia:*stai}d how an UBR 2'*pair is extraGte4.j v@ 
cimaidor a first ordar pair* Its and T ^2 paiaMsters frcaa 
«qn* (4*9b) are given as 



4. ' 

(1+as^) • 




(140S*“*) 


SubstH^ing m *- 1 ^ m get «»1) » +^1 and ** 

T^^C-i) » 0* Thisi^ at liss a 'otie* while T-jg has 

a zero, Ou%>nt is related to input and X^ (figure 4^7) 
as 

Y^ - T^^Cs) ♦ f^^iuy Xjg 


HftiMe, 


i^Cs) 


. a . T 

***€ • * *€4 


4^ 0 •» ± 1 


JSgps* (4«Sh) and (4,9h) give as 

■* i^mT^ y i ^Ca) 
1^^' #* vdiil ’•) 


(4*11) 


If h^(e> is %e> have a d^ris.-..!«e lover than l^C*) Wbu we 



1 ^'') 


+ 1) 


35 *« Hi 


0 < 0 < 1 


m 

(4*12) 


Hmm 8i5»r«cBipt (*) ^mims mnmUm w*r,t* ^ 


In slwrt, degree redt^tion of Itie transfer ftmotion 
by means of extraeticm of a firsb eafdigc* LHI tranftsissi<m 

*ero at E • H is equimlently ejKpressed as a *one* ireaoml 
from the transfer function* Blnoe, "to tmnsfer function 
is a structurally boundt4 fumtion^ sottetiffies a *<aie* can 
reeved only at mmm fraqiieiacy In this case ee 

extract a secoM order UH 2-pair described by eqn* (4*10)* 
Degree redisstion is achiemd if m let o as 


o 




1 ± 


n^\s f 


Q < Q < i 


(4*13) 


awe 


la eqns* (4*9) and (4^10), loeer signs are used if 
• -1* Also, if Jy(*) is unity at sT^ » ^1 then all 
.expressions can still be 'used aft«r r^laciag m by 

We now coiisider of the abowe synthttsis 


l«^iCMi!disre* 

WAMM 4** 


e«sider ttit following lossless transfer functim l^(*> 



H®r« If® find that H 2 C- 1 ) • +1* 

Ther®for©p ® first ®rd©r ISR J^'»^p®iir can he extracted.^ 
Uaking Idbe deriwiim^ m' get 



Q^m25 » 1. 


(1 ♦ 0*75 sH 


5 S~'' ♦ O.feas z~- 

+ 0*125 


nhich gii«s* 

1) m fl4/3 

Using eqn* (4*12) m ©btsin as 
Og • « 14/17 


Ii|(s), by ®qn* (4*11) is tl»n |iwn as 


*i|(*) 


*0625 


«*1 


whish is also a lossless fteictioa* as ®X|Nict®d« 

i 

Omd again* H^(-1) •• *1* 


Repeating the prooednre* m get « 15/31 
aoA Hq «• ♦I* 

The oasiMiled digital stmota^e is shoim in Fig* (4*8) * 



wg* 4*a 




KXAMPLE 4.4 


A »truoturaUy bounded transfer function Hj(z) has 

mity at - -1 and the function U given to 

te 

H3C2) • •», §^„ ^ 

^ 24 -• 5* " SJaT^ 4 IO2 * 


Suppose we first extract a 2 ikI order U® 2-pair at 
■ i ^/3. Eqn# (4*13) gives 

a m 1/4 

Also I ■ » Cos r/3 • 

H^(s) is obtained as, 

O 4. 

li|(s) • * ‘ "." ' >"" ' *"% s| " 

3 42* 

Kote that extreetion of a and order 2-pair reduoes the 
degree by two* 

Substituting • -dg we get iijC-1) • 41* «en©e eqns. 
(4*11)* (4,12) give 

0-1 » 1 / 2 , \ • 1 / 3 * 

fig* (4*9) show* tbi easoaded structure* 
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L I i- AIK 


-i' T- 1 . (/. 




i I. i (<f\ i.V K i 
I LK 

X. 

.... -■ ^ 


•> 



I 

i 




Fig* 4*9 Casca^M two^ipair structure of Example (4*4)* 


In Hm diaeossioxia mu IMi Jlwptir synl^sls* -we im^e 
aasimied the specif io structure deacritoed hy eqns* ( 449 ) 
and (4*10) is not sensitive to parameter o* Before m 
discuss its iii^leoentatioi% m introduce yet anol^r method 
of syniaxeslzing low seasitlirity digital filters* these are 
known as cascaded ortlwaonal digital filters* the name is 
so given because each tsiilding block is explicitly described 
by an orthogonal matrix* 

4*6 CABCAim REALISAtlOH W mmmmh OiaitAl FILtERS 

An alternative mth&d to realise a low sezisitivity 
digital filter is obtained by the direct factorisation of 
the lossless 'transfer Mtrix* lossy .fbnotion which is 
embedded in this lossless mtrix (secMcm ^^3} Is obtained 
by pfcqper terminations at the ports as shown in Fig* (4*5)# 

If the factorisatlcm is done into first or second ordin* , 
orthogonal matrices* nxm the strmtur# obtaii^d is of 
tim oasoaded orthogonal filter introduced by 0«l*^ttere&D«wildetg 3 

V« present here inily the basic ideas ©onceming ■ , 
these flltwfs# leaving out ^stalls ^ design* These will be 
iMlbd later while discussing tto mailed framework of all ; 
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low sensitivity digital filters, 

SteMUty of Orthogonal Fiitoff 

Ortl30g4^1 filters^ as oittlined in [26]# sr*® digilasl 
filters whose internal c<®^utational soheiae consists of an 
orthogonal matrix transfianBatioa, In other words, for the 
orthogonal filters, the realisatioB ma^ from lOa© inpat vector 
£ to the ©atpat vector j; is reppesentEid by an ortlxjgonal 
matrix. 

It is well known -amt an orthogonal transf<Hmtion 
Tg3t%Bmrm& the m>rm of izipat veotcsp. If is an ortivagonal 
transfmmatioii from 1, to j[, i,e« 

I - ws 

tlwn output noTB 

I*i - ([*] s) **(£*] S) 

• i. 

m (StBOa [T]'*f5j - I). 

An advantage of £t] heii^ an orthogonal matrix Is 
that it oan he represented hy very few parameters-, C©*»ider 
for iwm^le a 2x2 orthogonal matrix of the form 

rciwsi § ^1 


ii' 




flMB i^trix f(0) iti eiiii* inwlves only one 

paraiaeter 0*. Moa^eowr, if Is ol?se.rved tiiat if this paraaieter 
is pertnrhod^ my fro® 0 to 0 # o0p s*y,Xl the matrix TC0 +js^ 0) 
r<»aiiss orthogonal* In ©th«r ijords* no instability is 
inoi«Ted atltast by the pgqros^ter qnantiaation. 

Ihe^ above property is alsO' nsei- for a stable but 
very simple implei^ntatton of orthogonal transformations* 
fhe technique is knom as tte GOflDIC processing technlqtm and 
will be discussed later* But first m consider the low 
sensitivity aspects of or’toogonal filters in the context of 
l3*ansfer function reallsaMoa# 


Low Sensitivity of Orthogonal Filters 


Fco* lossless tpansfar. matrix it can be shown that 
aM have a coi^lementary relationship, i*e« 




" 21 ’ 


2 


(4*15) 


Differentiating eqn* (4,15), we get,. 

lT„(e^)| S| + lT2^(e3“)| Sj « O 


( 4 . 16 ) 

Suppose, that the tmiisfer functicsa H(z) is embedded as 
in the lossless matrix, fhen if H(s) is a selective 
function which is close to aero in the stopband and close %5 
imity In the passband (l#e, approximating Fig, (A,1)), then 
we have for tM passband, region, 
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- 1 

B®©®us® of ®qn* (4^15) m also tew 

T„(e3")| r 0 

Eqn. (4,16), tterofore, giws 

T2i(. 3“)| i 0 

This maos that m aybilamry first oriior parturtetion 
of any parasietsr 0 not only prasorws orthogonali'l^ and 
tenca stability lait also laaws tetevioor of HC@^) 
uteffootod, at Isaat in 1^ paasband region, 

(^r^bogimlity ^Uys an iaporbsnt role in tte desipi aM 
implementation of low sensiti^ty digital filters. In tte nert 
sectioii, m will ea^lain bow an ortbog^l matrix 'oan be 
lo^ed tgion as a eommon bnildiiig blade of all low sensitiTi%' 
straoteNis, We also diseuas tte implementation of ortte^osial 
natrioea by planar rotations nMch are oceipntationally stable. 



mnm fim 


im mmifinti smwimm t a yrnwim vmmoim amp 

iigj»!itjyAyios AspEc^m 

5,1 lifROP^IOH 

¥# imv9 »© far eoiyiidei^ mv&t&l iiethods of reaHslng 
low sonsitiwitf- <ligltal filters esoh of wMoli seems to stand 
apart frcmidlhers. These methods prodoee waipe digital filters, 
Gray rad Markel oasc^ded tSi t»«>-pairs and or'feo- 

gonal digilsl filters. The realisatioii scheiie of these methods 
differs so much, atleast apparentlyi^ from 0 m another that 
it is zsatural to interpret each of them as giving a imt low 
senaitivity stmotiipe, Tet all ©f these stnio-tares imvM a 
basic feature in ocwraon ]Mimely»1hat they consist of sisiple 
building bloolES as pmto^p^ l^mtt are commoted in a 
ceiscaded fashion, MapeotWf a iwtiiw ©f l©8s3es®»ss is 
associated with ©ach ©f #wse str«w1»«6. In’ the case of 
lE^s ai^ Gray and Marteel* s lattice filters* this losslessness 
is described in terns of pseudopower which is related to Ihe 
incident and reflected signals* twO'*paii^* the transf^ 

raatrlK in which the trai»fer fwa©ti«»i is «Bbedded is lossless, , 
Finally* for the orthogonal digital filters* each building 
blodi iniOTospoods to sa orthogonal transfcrmatl?Ma which preserves 
the mm of input yector| . 





13 


In vi«w of the ahoire ai»ilarittea,an atteiq^t is maae in 
[ 0 j to find a link betw^n wrions realisations mentioned 
earli®*** section# ^ elaborate on this laiified 

theoretical frwsework within which a close resemblance among 
these structures is clearly eslnblished. Later in the 
dis 0 ussionS|, we consider certain issues connected with the 
implementation of low sensitivity digital filters, Ve discuss 
a shahle tout simple sches®® known as CiSiDIG 

techniqua [273# that has toeen Introduced for the orthogonal 
l.iBplementation,We also touch some fundamental concepts 
meded to understand VLSI iBg»lementation of these filters, A 
geiieral ssrnthesis proeednpe for Vlfil design of these filters 
baa.d on «» «ort prowntwl to [ 7 1 ha, al.o bean toel«lad 
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(1 . o.-') • *22 (V ♦ 


■12 


jT - au 4 ^>. 


(5*1) 


To roolise this 2-f«ir ag aa orl^ici^oaal matrix^ im 
iSx»t extract a delay eleaeat fro® it as sbcim in Fig* (5*1) * 


X. ■■■■—^ 


L j 

i - PAtf ?. i 




«. -..I I 




A — ■ V 

IT ^ 


Xt 


Fll# 8#1 Kxtmotioa. of delsy-^leoent from an UK SHpair, 


figure (6*1) ean be locked ii>cm as a digital 3-*i»ir 
uliose port 3 is terminated in a ^lay element,^ i*e# 

Xj m s"**^ (5*2) 

Siig^se that •mm digital 3-port lias tl® following 


r <if 




^((»e 

*3 



i 

m ^ 



whcr®, 
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ifll 


Mi 

^21 


Mi 

^12 


®22 


• 3 




31 


als0« 


£ - [*1X2]^ T - [Y, y^] 


5uto«ti”teitliig ©qn* (5*2) 111 (5«3)t we ofetalu 


Y* • 


T "^X 

m » 


^ (1 

f h© transfer matrix Y(e) is gi’wen as 

..0^ 

Tim) m [T ] 4 t f ^ 

** a^ (i *# f* E^) 

% 53 * 


(5*%) 


(5.5) 


Equating term I»y tsnsf C5«1) and (5*5) , im list 

folloicix^ fora &i ■ttsa transfer i«trix f* of' digi-^i^ 3‘-^irp 


T* 




o 




# 


o 


(5,6) 


that if is scaied by a 
matrix f(E) 


It is obsarvod from m 
faotor « and a faotsr 1/« 
is not affoctsd, ■ 

If m tboosa « *• J$f sqn, (5*6) is 

rowpit^ii as ■ . 



T6 


I**® ' 0 

0 0 

On® ©an verify tiiat tills matrix is wthogaisal, fSias, 
a digital UBR 2*iE^ir be realise4 an ortix^osal matrix 
if the delay is extraeted-firaii it,, AalRial imiilea^taticm of 
eqn» (5*7) will be discussed latert First m consider ttm 
case of nave digital filttfs* 



(5*7) 

0 





In lAve digital filters ^ adapters are the btiildii^ 
blc€^ of ladder stmcteei as shown in Fig« C3«©)* For 
example, a series adapts C%pe ®2) is given in Fig* (5*2) # 
Pwt 3* which is terminated In a delay, corresponds to a 
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> i 

r 

«» 


»i ! i 

1 


■^1 

i 'll 

f 

Bg j ■ 

' i 

o 

T 

-1 j 

1 

1 

1^2 

1 

B , ' 

. 3 ^ 1 


1 

'* A , 

f ' 5 ; 


If port 2 nave quantities, ^ mad are scaled by a faot<»? 
©f V as shonn in Fig* (5*5) and **«;*'^ is extracted fro® 

«4 

port 2 Instead of z i*®* 



Ai 

Bx 


> 




&FLAY 

eVTRACTtCl 
St Bies 
adaptor 


-W, 


B 


mLttW 


«:i 


. I. 

Aq 


fig# 5#3 Modified, struetw® of series adaptor of 
Eaonple ^ C5«2) « 

then ifo have transfer isatrix In eqa* (5*8) as. 





%r^i 

§ 

\l 
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Bqn, C5#9) Is same as eqa* (5*6) with » o* Bence it can 
also be implemnted as an cartlKjgoiml matrix as in eqn* (5*7) 
aftar proper scaling, 

Om mn similarly sh&w that a parallel adaptor too can 
be realised as an orthogonal iiatrlx* 

5#5 d^HQQOMAL OF OSGAIED UTOGE fII3!ERS 


We now consider the case ol Cray and Marlcel filters 
whicli consist of a oasoade of SHport lattice or ladder 
struotures* Chain matrix of this SHport is given as 


mi 


zB, 


hrt-1 




■ 

A 1 


m 


■ B ! 


' » 1 


(5*10) 


Since the transfer fiMctKm is realised by nsing ratio of 
and A (a » 1**, ♦*!?), chain parai^ters can be scaled by «iz. 


as shoim In Fig* (5*4) « 


A.m, 


:Sni*i 






'A 


Ti-4<4^^6 


m 


■fi . 


Fig. S> 

figisre 


strt^ture of 


Thi foaled chain pa»aiiet»« «?® 
Jk. e 0 • B • C 


C5*11> 
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Eqli, (5*11) gives ttm trsimfer Eoatrtx i^raseters 


ss 


T mm ^ 

ill4l * JK,^ 


r 


‘11 ■ T22 " *12 • *21 - / 1 <H^ 

TIttM ■m G&n write* 


Vi 





B 

m 1 

*« 

m 

J 

a 




(Bin) 


(2^*12) is realised as ia Figi^ C 5 * 5 ) and is known as the 
.j|Q 3 aii,liJMwi. Pray ...ani. .ffail Ettl, Xadd^ heeanse it is an orthegenal 
«tri 3 t^ 


MiWl 


c5> “V- 


Am 


?e 


t 

4 

"tCI" 


VPk 




ma >!^ T '^^7 Bm 

l<rr: 


Fig* 5«S lemaHsed Oraf and Mfdiel ladder structure* 


Disoussions in m&tim C5«3!} * iB»B) suggest that all of 
tbs oonsidered low sensltiiritF digital filter stru®t*res oan he 
realised as a cascade of hlo<^s whose transfer aatrlces are ' ' 
ortliog«mal* It is the cridiofaaalltf' of Ibese aatrices that 
a^es ns conclude the fol'lowingt Afferent .aethods o^ 

. ,se»siti M1lf-.dlJE itil.,lilt^ give different struct 
.■ttoiy...h.saa» tia 41 y l>e.long ft f|3.thrff« 
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In Wm next snb-eecticai w XmM at Bm& of the problen 
iBjpleffientatlon of low senelttislty strtiotnres [23 ]* 
Comnonttonallf^^ they isan be iapl^i^nte4 on any general 
piirpoi^ oo^nter or a de^icatei digital aignal prooessoTi, 
These ii^leiaentation® |n*ovide flexibility in the design* 
litardiiax^ in^lenantation is inref erred when speed is 

Of prill® oonsiderationf Slmilarlyg when mass production is 
required we med to design low sensitiwiiy digital filters 
suitable for VLSI iip^lenentations,! Before considering the 
VtSi design requirenent^g m shall see horn an orthogonal 
Mttri% the * oore* in all stmoture can be Inplenented by 
successive * planar rotations* « 

5*6- ORTHOGOHAI. MATRIX Ml FUHAII ROTATIOIiS 


As it was said earlieri^ an ^t^iogonal transfoxmtlon 
does not change the nor® of input -yector after transforaation# 


It win be shown now tdiat it nei^ly “rotates* the coordinate 
axes i*e« tlw output vector coordinates can be obtaii^ by a 
sii^le rotation of the plane lav J* 

Consider for exasple equt (5#12)* If we substitute 
m cos 0p then the transfer mtrix denoted by T(0) is of 
the fora 



sin 0 

0 





m 


Ae sliotm in Fig,» mtarlx earn Ise 

as a ratstitni of tkm pleuoa forced ^ mmm asii 


; . , 

"■=».v, 

Figt 1*6 @f cr-tasog^l aalrix by plamr 

ra'ia-y.m# 

Sindlarly by siibstAtutiag # C f«r* W£0is) m m& 0 in 
©qB* CS*7)t «« ©btain trsaefar aatrix f(0) as 


i sin^ 0 e©s 0 

I eos 0 0 

j^slB 0 0 sin 0 


sin 0 ©as 0 
-*«lll 0 
-NES©S® 0 


CS.14) 


This ortiiogonal loalriie ©an be faetisrized into a |>ro4iiot 
of twa siatplez* ortbogonal wtri©eS|f eeeh of tfblob «5orres|>o®^ 
to a planar Fotation similar to tb© one shoiei In Figj> (SaS) * 


The faotiarisatioii is as folloiw 


7(0) «» t2 (0) f^(0) 


l oos 0 

sin 0 0 


0 

1 

6 

0 

0 1 


sin 0 

0 

oos 0 

sin 0 

*«OOB0 " 0 


oos 0 

0 

*i^fn 0 


Ci#ib) 


’ Fig* (5#7) aboi® b©if fee .I«lEe4 aa 

a ^abinatlon of two planar rotatiana* ■ 




Fig* Stf Siae€s«s»iip« plBjasr rotaticnaa for a general 

siatrix* 

In (4#6 ) it was shown that the orthogonality 

ordk^ 

of an ertlwsgoml matrix is not 4is1whed hy any firsVpertnriaa- 
tloai of the parameter 0^ This results in ti^ stable config|;i* 
ration of orthogonal implementatloii of filters* ve xmv oonsi^r 
an algorithm known as CORDIC alg€a*itlas whioh is very suitable 
for implementation of orthogonal matrix by successive rotations* 

5*7 THE * CORDIC* AhSQRITIM 


The primary purpose of IMs algorithm £27 ] is to oim^tite 
the output yeetor B which is obtained by an orthogonal tran®** 


formati^, TC®)* of the ii^ut wctor A* This is given by 


X 

n 

lii 

i(«)] 

© 

n 

m >d 






# 

' * 


t*n 4* ■ t*o 



&m ® 






Tim ©©aputatifim f«r Itiis is ©ajrrled out by 

smmB 0 i'm plmmr rstaiieBas of tte iisput wotor* Tl« sallsal: 
festtirt of tbo alg^ri-ISb® io that only sMfts and •eross- 
additiiPiMt* ajfs roqnired* TIsis i^swlta in iiie aiDaerleal 
ataMIity of the ©cwg^^ntatioiial sohM®* 

The algoritbB 1» as folloifS, Si^jpose, that o is the 
angle by whioh input iwoto®* ^ is ^ be related to get the 
oulp'iit 'wcto®* let a be est^i^fessed as a lin^r ctmbinatiO'ii 

of 'aiig.les ®j|f I 

H 

am S X, a. (5*16) 

iiO ^ * 

Here can be +1^0for i » 1 to n and they are also called 
Idle Arc tangent radices - (ATtl) digits* 

If n is the length of the B»B«»ry register and each bit 
denotes an angle a^p such that 

O - Tan^(2‘^^"’h, 1 > 0 {5.17) 

i ' 

then each nonzero bit denotes an angle aj^# 

Also* » 90®* 

Consider mm the sittation sheen in Fig* (5*8) ehere a 
yect^ OAj^ is rotated by an angle ♦ 
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fig# i#B CCBDIC appr^xiBjation of planar rotation. 


If angle a, ia small enongli then the dotte<J arc can be 

m* 

approximated by a straight lin®. In such a case, coordinates 
of are given as 


* 2 - 0 -'*) 


?ii2 




(5,18) 


Eqn, (5 #18) can be very easily implemented by nsing simple 
shift and cross add operations# fbe typical configuration for 
this purpose is shown in Fig# (S#9)* 



Fig# 5 #9 General structure of GCBDIC iE^lementation, 

It is noted in Fig# (5,8) that any rotation tbs 
Magnitude of vector OA.^ is increased by the factor 1+2""(i”‘^)2 



SjLxim Imngth n &f tlie' re, ister in flaetd ttee ewimll 
lnor«B@»ta3. eirer k t» & eenstent* fM» is given m 

k * S (5*19) 

i«»1 

To illtistmter the nlgoritlm m eomlder following 

EXASIPiai s*i 

Ciii^iite coordina'fc&s of P* oktninod afier rotiaiing P ky 
If® oomtorolodcwiso naing ftaRJDTC alg«irita« for 8 kit iiapiefflantatio 

Her©! lasing oq»» (5*17) a^*s are given a» 


% 

m 

90^1 m 

m. 

14 

,036® 

“4 

m 

7*125®^ Og 

* 3*5763®, 

®6 

m 

1,290® 

a_ 

7 

w 

0 * 8 ^®*’ ttg 

«K 0 *^^ 8 ® * 




For 

“ m 

« 19 ® we mm i»ve 




^2 

w$ 

<*1i, 3v^ w 

— 1|, Xq « 




0 

m 

^<1 * ^3 *• 

Xg * kg 


kj 

w 0, whiok give® 


1@*9f!lS tlie ecspntation is oanriei out as foHoiii t 


0. m 
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1 

nmiim 

h 

Tj 

i 

© 

0 

jmmm 

,01100000 




♦*00110000 

2 

♦1 

,00110030 

,10010000 



♦.^EWIOOIO 

•*00000110 

4 


,01000010 

,10001010 



*,IX>000001 

•,000CK5000 



,01000011 

,10001010 

fhiis 

the conpu'ted result is C*26l7|i 

#5390) where 


actual result Is (#239^^ 


¥e ha-ve seen that it is jmssible to ii^leiieiit sm 
orttiogonal B«trix by a irery stable planar rotatitm B>etlio4*^ 

This method is also suitable for VLSI impIeKUtaticm cm account 
of its •modularity* ♦ In the following disoussions we look at 
st»e of the requlrea^nts that need to be met for the MSt 
impl^entation of these filters* 

6^ RBQUIlSMSIiiTS FOR VLSI DBSIGi ' ■ . 

With the recent advances in mier^leotronios* certain 
digital filter structures that can be realised on a single c^p 



hav<e b®ec«e more attraotive tliaa otiiera* Certain terma need 
■to toe defimd befcap^^an afqpreisiato ■tte mm f^uii^ments for 
fliSl daplementatlim of toeee filtora* 




tim most important oonsideratlon for a stinictiire to toe 
VtoSl iaq;»lem®ntatoiOi Is toat of aodnlarltv * In other i«ords^ 
toe entire stniotiire must toe a i^gniar interconneotion of 
similar p»E»i^8sixig elements,^ it is. also, oalled an array type 
aroMtoo^sPe and toe easily fatorioatod on a single otoip* 



ConneotivitY (Sys tolio liature) 


One prefers in VLSI design ttiose strnetures wtoicto oan 
be dram in a single plane vitto no *orosswiPes* • If in the 
entire array^ each processor toas only -toe nearest neightoonr 
links, toe stmcture is oalled systolic in nature as stoom in 

Fig. (5.10). 


■f 


j j-vj i-i-j ■ j->; 1 - j 

1 ^-V-| ”* i 


I f'i p-'i |- 

1 U-| , 1 

---j 


Tivo DfMC N LION C SfiSU AR C 

A ONI- DtiMt N LION' OY'Tour. At,'';7-\Y 'YlTolic /£k?r^Y 

Fig«j^0 'Exai^les of systolic strnotaiires* 
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aiT«y type stmoture^ the proeessiag ti»e for each ii^ut 
data depazida oa laagtfe of array awa toagh it la syatoHo 
ia aatiap®*, This is beoatstae of Idte fact tl»t the ii^put to 
procNissiiig aXesaaiit is awilahXi only after . pr-ooessiiig by 
its eXeneatCs) is oo^lsted*. Zmh sttm^astee are 

^ .titg.liiieabl» digital filter is one ia which the 
mm$mm throog^nt rata does ^t m the <^d@r 

of filter which is dlreotly proportioiial to length of the 

ar^sty# 

It is also desirable to design VliSI struotnres that @an 
iiapleiaent digital fillers with a ^wlety of -te^nsfer fimotiom 
Of different orders ^rely by ©hanging certain paraaetars* 

Such strtictaires are called nnifona in nature and can be i®ed 
to process several input seqmnoes for different transf^ 
functions* This resiilts in increased throughput of the 
structure and increased ’active* ti«s of each processing 
elment*' 

Having introduced seme is^ortant co»epts of VIST 
iuEpl«a#ntationj> w# consider the suitability ©f av&ilable low 
sensitivity structures for the VlSl design viewpoint*i[ 2® ], 

53 yisi coMsimATKii foa low ^i^iTmTY MQitkh m.wsm 

All structures of %mi sensitltitf digital filters 
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TMs iBcXutos mm liigital filters als© ewn 13i«tigb they are 
rea3Li8e4 hy an Indiyrect aethed^ Sys’h>lie mture to© is i«reseiit 
In ®a©h of these fll-^r stiwfaihs,, 

PipelimebiHty, l^me-ver^ neeils to be investigated 
fiirHajeri, For ■ 1HFS| in whiidb f#fcllel series adaptors are 
alternately oonneoted (Fig* 3iMll§ o«ti>«ts of any adaptor is 
ooi^nted froa its ii^wts* These iiiJiits intum depend on cwcrest 
of ad^aoent adaptcn^s* As a resulti^ the struotire is 
non**pipelineahlQ* mwm^ the larger is order of filter 
transfer fimotion* is the prooessiiig tine for an Input 
data* On the o'^ier handi Gray and Mailcel and 2-^lr 
structures are pipelineahle* This is so hecatme of a ^lay 
introduced between successive bloc^* Thus processing tine for 
an input sample is the same as that of one processing el«wint%, 
Such structures are also called svnchroimus because each element 
works in synchroniast with otJ^rs* 

None of available two sensitivity structures is^ 
however^ uniform in mture* This is so because given a transfer 
function* by merely reodering the cascaded blodfes one can not 
achieve different realisations of the design with different 
parameters* In the next sub-aectlon* m pi^sent an orthogonal 
digital filter structure which is both pipelix^eable and imifcni* 



In a recant worit {^9 3 a scheme^ to be presented lare* 
has been discussed which has all Wm featweas'iif fI3I 



r^quirwieiitsi, In -^lis design ■{^he tesl© l>nilding bloc^ hns & 
^0 tmnafer matrix and Is sh<n«n In Fig^ (5*11)# 
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Fig* S*11 Structure of an wthogonal digital filter* 

Fm* &mh filter tte transfer siatrix T(z) is (n'tbogcnsal 
and the input output relationshiF is giimn by 

* tCz). 

C5,^) 

PP* ♦ RR* - QQ* (5«al) 

Inhere %(«) *• iCl/z)* 

fioe, for a giwn t»n»fer function B{z)^ thaire are 
seireral choi^ts of T(s)* St;ppose im wish to use 


yj^Cz) 


: Hz) ] 




yt3(z>| 

i 

1 

: J 

RC») 


the orthogonality of tCe) giires 




C5#2a) 


3K»5 • f 

Snell tiaat Miz) is ^ transfer fimotiozi £f£« ii(s) to yCz)# In 
tills tho followitig relatidnsMp Isolds - 

p(s) • Miz) 

«C») +: k^iz) m niz) (5^25) 

Wmpz ^{z)^ B{z) ars tl»i niaieriitor as^l dencsedBator polync^ials 
i^j^etlvely of transfer fumsttoa E(is)* Eciiiaticni (5«23) snggssts 
that by serely ehstngii^ pararoters Ic^ and diffearant timnsfer 
functions can be obtained, Furtbar# tl» ss®« transfer functliai 
can also be rgaUsed in acre than oi» way* f his result® in the 
'uniform* nature of the filter stnwstur# as described earlier* 

For illustration* we consider an ernsple* 

Eealise a transfer function H(a£) « p(s)/e( 2 ) for the followi 

cases 

(i) • 1* 11:2 aiid(ii)ki| • k^ * ^ 

£feat-l. 

SubstituUng • 1 and k^ • 0 in eqn* (5*23), we get 
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P(*) • p(z) 
qC*) » #(®) 

If we repiece R in l»y f tl»n we mn writ# 

. ee» « fp* 4 f % 

wtsicfe i» saa* as the feldkeller tm* (1#2) few whioU 
Stantoi solutions are awailafel© if H(z) is a gemrallse4 
transfer fnnotion (A ^1)^ 

(ii) Case II 

Snbstitnting ic-, • » 1 in ©tn* (5*23), one obtains 

P(a) » p(z) 

Q(z) + R(2) « ©(*) C&#24) 

If we ehoose QCz) • (eCz) + fCz))/2, RCz) » («(&) 
then substituting in eqn* (5,2^) we get 

efi^ + fa, •* 2ppj|. 

IwDleaentaUon of Transfer Matr^y 

Qmm Hz) is obtained using eqns* (5*20 - 5,23), it fs 
possible to faetorize it by an algoritha iKi3 ) into sl^l® 
orthogonal iwiilding blo^s [f. J as sJ»wn in fig, fiie 

chain matrix of each block can be written as a product 
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\ “ ®i2 ^ ®i1t ( j » *ar« girm fey 
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1 <m|£. 
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© 
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0 ^ ’!'**%. 


II' 


2 
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ia 


— %. 


1 0 
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"12 
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^’Tlk 


12 
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C5#25) 


Vbere kj^^ig (5 •• tfS) ai^ ©btaineil as In algorl-toi (A,13)|j 
Fig* (5*12) shows how aaoh ©an be idealised as a easca^ af 
two normlised Giay aM Markel's ladder stmotiares of Flg^ (5*5)* 
The pipelineability of the structures follows fro® the faet 
a delay element Is present between any two fenillding feioeks 
(Figure (5*11))* Equations (5*25) ©an fee isq^lemented fey 
planar rotations* m shall pre^ntly see, a diff^pent imm. #f 
CCMJIC algcrithm iapleaentatioa results in ti» pr@iw>##iiig ®f 
parallel sequences in this stniot^re* 


Fai^llel... Fr ocessing ..of .....IjiiPiit .ieduei^i 

As stated earlier in the CiWIC algoritlai an gula r . 
relation is aohiewed through successlire steps* say n In maabar* 
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V ^ I 


'■> 'j it 

4 ^ 4 Ji k, 

^>'A 
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l‘■,: 
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‘■-A 


/■ 


I /i' k, 

©— <1 





l3Wt»ad of th« ooaflgxarmMon of Fig#. (S.# $ )# if w© ti#® 

•forward going* sbifts and cross-adders a# sliowii in Fig# (S#15)# 
wa !»£«*•■ that aaoh '* shift-adder* Is sstiire only for a fracti. 0 n 
(1/m) of the emtinr |)ro€»8simg period f # ifencs if tmf fer 
registers of length' m are used for 'parallel input sequencses and 
AfR constants for differsnt transfer fuMtions are atcn*ed in a 
circular registers of lenglti ,% ttjen it is possible to process 
several input sequences in parallel.* 

In the discussions* m have ®imm how different 
structures of low sensitivity di.gital filters are realissthle 
using orthogonal building blocks* ■ fhese blocks are iaplsasented 
in a stable manner by successive planar rotations such that it . 
is least sen,sitive' to the perlwhaldons in the parameter values# 
In a broader perspectivei ttieory of low sensitivity digital 
filters is also able to offer cfi«#etent designs for liii flSl 
i^lementation of digital filters which is not guaranteed by ' ■ 
the oonvemtional design methods.# 





tn ttiis mxii. we imm ^smjyi«e4 tim design ef lew 
mmiU.vity digital filtm using different Isinds of 
Siisulatioii results B) for -^ee 

etnietwes ham been shown to destcmstrete their sigierior 
per€Orneut»>t ower c«aawenti<»iel settK^ds of iog^leeenting IPt 
digitel flltiors# It fans also been sbown that a general 
Iteorr based on the oimcept of ortbogojmllty' provides a sore 
eonorete and satisfao^r^ erplsnation for the working of tbs 


'varioim low sensltivitir s'^raotin^Si 


Imw obsemred tbat 


many of these structures are also suitable ^r VlaSl |japl«»in* i 
tation* Based on tbei^ facts » it is reasonable to conclude 

i 

that low , sensitivity li^lementaticffls offer a aore attwictive 
alternative conventional schei^s in asst hardware appliea** 
tions* 

It was originally envisaged that ttae design ©f FIR digital 
filters will also be investigated* for these filtscs^ tbs 
realisation nap frcm ir^ut sequenoe to ou%)itt se<|ueme can 
be considered as a linear transfonaatlon ©heumcterized by a 
cyclic transfer matrix* It would» be Interestii^ 

t© see if it is possible to ii^se certain conditions on the 
matrix, say orthogonality, fcr a low sensitivity iB^lementatioo : 
of ' bloc^' ii^ut sequences* One may possibly te able t© 

©btaiii vector reallsatiiisi stfucttres analogous to the parallel 







APPgHPn A 


A.1 QKW^AllBm FllffiH fwc^om 

Generalized flltel* functiozis are used to approximate 
14ie i<2Mil ireqwnoy respcKOsen, A® Meal loe paae fr®<|ii©zioy 
ra^xwiMi ia ahom in Fign 

ImCjiiOI 

J-0 \— i 


__ % CW') 

Fig« A'«1 to ideal Imp&MS irequeney re^^aaei, 

The transmittan©® of a filter is expressed in terns of 
the generalized filter functioii K($^) m 

1 



For n order Butterworth filterp 

(A,1) 

%diere n i order of the filter^ G i a ©cnstant 
or* 

I “■’ ■ ‘ ■ ’ 

jf k and are maximiza and adniBRai pusslmM and stopliaM 
at-ientiations respeotiirely then* ■ ' ■ 


a ^ 






(*, 5 ) 


iflioini 




m ^1 I K ^ 


1©' *w'f 


CA#%a,^) * 


yalog l«a«t v&lm ©f a and meetUij «» »t©pbaM f^nirewant 


d 




I© 3^# 


.» 


(A#4©) 


f«r ©C») giirea, {^©Bsid^lag ©niy 


zmrma for ataMlity) 


©C«) 


n 


°kM * •!£ * «P(JC 

fe • 1 




Ctxebvahev Fll-teri^ 
.tki 




n 


{A,S) 


For ordor Chebyshov filter. 

|K(3U)I2 c 2 ^(q) 

Idler® Tjm) is tbe n’’^ Chebyahev fraction 

t^o • If ^*1 • ^ I ^ ^ 


(A# 6) 


Cosb^^Cl/iC^) 

n > — ""i_j|i I III ^i.„ 

Cosh“^(1/K) 


(A*T) 


idiere K and C are gi-wen by CA,4a^b,c) abow 


*(s) - * 0^^ (■"•fc), », 

• Sinh (1 SJjflr' 


♦ 


m} 


w 


\ m Qmh (i 1) k) 

fw Ic • 1 ‘ to ii an a^roijrlate oonatantiv 

A«2 im:i 


A p^ynamlA}. witii Ita roots strictly In tiie IBP is a 
f^ysMdLal* ' if in addition l^ing hunrits it is 
to tow siiB^la raroas cm Wm ^ asls tben it is 
iuMpifits, A f o.8itiw r«al tmmUm is cam imich has 
foUoving frope^tiss 



(a) It has rational ooaffioionts 


(h) Deiamizmt»>r is hnrwits me* fflixiife4 hnrucits 

(c) Real part of functicm (wi^ s * is non-^gatiw 

(d) Residues of the S^^-mdLn poles are real and positive. 


EXAMPLS A.1 

Cheese if Z(s) «* is a 

S' + 2 s 

(a) is obvious# s *• **2 is in liiP and s •» 0 is a sia^le roo% 
cij axis# Thus denomlnalxn* is modified hurwrite and (h) is 

«tlsfied. Also. 

RgCZCd^)) • " y -— g I which is non negativei being even 
function wilii positive coefficients* Hence, (c) is satisfiad.^ 


Mm sZ(«) • 2 > 0, satisfying (4) 
s-^o 


Residue at s » 0 *» 


ii 


h»m« zCs) is a FRF* 

la as'^fCHik iSmtayp flFs “tilMEtir iiqpsrtsz^Qfs frm 

ttm psstilt tiist a tiaistloa im pmittm raal itt it is r^iisalils 
as aa lanittaacMi of an RLQ! aetwca^* 



Im #iis am isnit^i^ fmotiom is rsalisad hy 

pole remo^Is# Firsts all poles om ;p axis 
laoliadiag t^ose at s » 0 amd s m m are 3?aiiBowd« iMs is 
giwm toy equation (A, 9) and glws Foster* s syn-&esls if 
FjCs) is zero, 

Hn) - -f ♦ ^Tr * Fj(.) {*.9) 

i 

where are residues* Ihe network is realised as shw^i 
in Fig, (A,2), 



Fig* A*2 A Foster* s »etw<8Plc* 

FjCs) is ©ailed the miniimm funotloB ainoe it oam ix>t toe 
radmoad further toy extracting ei-tt^r an iixluctanc^ or a ©epa®it»iwe^ 
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It is i%iil48«4 ty rwsvini a ©©ustRst m that iiie ri^iaiag 
p«?t is still a IMS is. illMstea-^i fey aa axaag»li.. 


Ss # 1 


at s « sa 


s— >a©’. 


...±.....fm..±. 


s(6s ♦ *1 if 


if tar af this pale w» }mm^ 

ZjCs) - ZU) - 3» - |f-|- 


Siaee 2jC©) m 2 ^4 2^ 

Ej(co) . This leaves 


Cel m ne resaye a minimum eaastEtut 


2 j C s) 


4 a ♦ 

iSs"+'' 



Henee ZCs) is realised as is Figt C-ifS) 

— A/vvw — 1 




Fig* ^#3 Fester’s aeta^et ®f 


(■t #2) ft 


la the Gamer* a realisatiom, ^ W , is redueed by eae 



dcjr@«| tioe a pole at # • co As roiowA ti'Bier JCro* Z(») 
or trm Y(a)* It ganarataa a ladijtof as m&m i© Fig* (A#4)* 


L 3 





Fii* A*% A Caner* s laiilir mimsk^ 


Q&m&idor ZCs) *» 


Removing a pole at s •» oo 

z(8) .8 4. . toaLt-as* . B + Zi(») 

s + * 

Removing a pole at a • oo from Y^(s) 

Y«(s) » T^jCs) - ■* 


Repeating this procedure ne get tiMt entire networic as sitioen in 
F ig 0 (a *5) * 
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FAg* A*5 Ladder aetworlt of Example (A *3)* 



for a spaaatria lattljee oetwaile 


^11 •“ ^22 • 


* %1 ** **" 

imtmixmVLag mai^tim& for w E 


♦ I^E - 


Vj + 1^ 


0 


We oto-faaiii 


(Zb-^a) 


(2R + Z^+Z,,) 


(A.IO) 


and 


a 


2(2R*Z,*Zb) ‘‘ 


(A.t1) 


Subatltuting eqns* (A *10) and (A #11) in taie expression for 
t^s) vm get . 

CE't’Zji) (E't'Z^) 
m 

'2 


u(a) ■ 1/t{s) - I 


Slnea Z^ and Zj, are losalasa and henee odd functions, «a taava 

,(R+Z,)(R+Z^),,{R<«J(R-Zl,)^ ^ 


ii 


iCi) 


*7 \St 

iiii |ji !ii ) iiii i i I I I ii jK i ip ii lBlMPa 


X(«)E(-c) 

(rR » Z,/^) 

R(Z^-Z,) 




I# 


"ysg'" 



R^ - Va 

R(Zb • Z^) 


ifliiiiliL giim the ex^essions fcr aM Z^ in tezies of V % aM 

A,5 SmSILIfY of MICROWftm 

I 

For s » 


f « Tahh(^t/2) 


« 4fan(wT/2) * 0 



Thiis the entire ijj-plMm is mppeA onto a strip frm Mi/t 
to 4-n/f In "toe e-*plane« Horeowr, 




^ < 0 < 0 < 0 

t ■ 

ttimtpr#, a pole in IHP &£ ■ plmm to msxppe^ im ^ IM? atrif 
in mplMm presarviag liia ataMU% : ^ . aoaditian 

ttMi 



If Z(«) is a MF and lCs)/Z(1) is not aqual ^ m &p 
i/% itetn a unit elsiaant of iipedaiiuQe Z(1) can alnajs be 
aartsraoied trem zCs)^ laairing a 2 ^( 0 )* If Z(1) • Z(-»^)^, 
then Z^(s) is lo%^r by one degree* Z^(s) is glmn by 


Z^(s) 


2(1) 



sZCs)-«Cl) 


(A^12) 


Tim theorem holds for adaittanoe f^naction TCs) as mJX 


with Z replaced by Y* 


A,7 lEUY FREE U?OPS 


Delay free directed loops can be generated idiile 
Intercoimectlng an adaptm* pisrt to an element cr a peart df 
another adap-faMr, Both capacltanoe and ix^nctanoe involve a delay 
element* Similarly a resistance is as open clrcnit* fhiJ«t delay 
free loops are never generated idienever any port is connected, to 
RDC elements, Howeveri^ they may be present when a port is 
te rmina ted by a ’short* or ’open* cirenit or connected to 
another adaptor as shomi in Fig^* 




trm Ic^i* 

' . fo ainodd '^sa ialaif £k*@e m mn mTx.e om i>f 

l>©rt* to be eounected m mfht&tXm trm^ Tims t&t aja 
n«p<wrt It' m mhmm 



Rjj “ ♦ Bj ♦ Rj^ 

(**13a) 

or. 

°n ” ♦ <*2 + .... On-1 

(*.12a>) 


Here eqns# and (A^ISb) bold for* aeries and 

parallel adaptcjra respeotlwly* F<^ tt»ae cases, la 
made equal to unity. Port n thus beoMes refleetion firee* 



For a ladder each adapter • ahacrha* &m "reaottue 
element only (i#e, L or C) , third port resistance 

©an be arbitrary ahosen so as to «k® it if reflection free* 
Third port of ttm last adapter is terminatid in a resistaneei^ 
beiMS® it is an open ©iroult* This gites ef^utability of 
any ladder waw digital filter. 
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A Q 9iP*&*P*PR.OTWl MAWTV WTSU UtW 

»0y ..HftHuA, .fyjf ...mJEMM.. .mW 

for o^araetrio lottioe of fig« voltog« 

and currant vactora are related by 


■■ 

1 

« i 


2b - 2. 


2l 

*a 

• a 

r 

.... ■ 

2«* *b 


1 

w 

L_ 


• I - Wi 

'V " ...I , ' . ■ ■ ’ 

Using the mm pareusetersy m urite 
i • X ♦ R[U] i 

f » X -rM g 

Vimr@, A m [A^ f • m m R 

[U] denotes tbe identity mettrix* X i ®**® related to A 
and B by 

V • 1/2[A + 6], I m 1/2l[A ^ S3 CA*15) 

Substituting eqn, (A * 15 ) in (A,i4)aiid rearranging tiie teras, m get 

a - -(W + [zJ/R)"’ (M - [z]/R) t 
- [s] a 

EjOpanliDg tZ] and alnplltylagt Zlni that 

aj, - (Zj, - Z,) R/(Z, + R) (Z^ ♦ S) (A ,16) 


If ar® aafinea m ^ " 

***" **' ** ♦ R) 

fhm it is ^sii|r aiiawii that 

GRAY itg) JIARKEI»*S iimCRIfm 

F<»r H(a) as la equ* (4*1)^ Hie aigisPiH® is as felioifa* 


t*1 

l<et a 1 ^ « n ’ 


ta 

I * A^(i/a)a"^* i 

A„<e) Witten nlth rewrwd 

of coeffictenHi 

1,3 

Vi » “ V" * 

last coefficient ©f A^Cz) 

1,4 

\ * " Pb,« * 

Last coefficient of Fj^Cz) 

1,5 

Aa_l(z) ! - LV*> - Vi 


1,6 

8 = Pjz) - ZB^(b) r. 

1.7 

m 4 » ffi -* 1 


1,8 

Go to Step 1,1 if m > •» 0 



•» 1 and Pq • 1 Sf® defined. Because of Step 1,5 
is unity for all », Usiiog recursioa m S1»p 1,6 aii!l adding all 
terms m get oqa, (4,2), 


THE umcE km uxsm mmmmm 

frm step 1,2 (A#1'0)if‘i« get 


Vv*) - b^i(v«)]/<i-h^) 

ZB^,(*) - \M {l-»2) ♦ 

$ttteM1»tlll% 9mi <A*1@) $M CA#I7) wt muin m 

Vl^*^ *• ' * \ . CA#1f) 

t»m iifkift vwAng CA^Id) ' tmi iA^i9) m g«t 

\ iA^O) 

i^e., CA#19) (A^O) gif® ids© lattl®« stnaetiare relationi, 

A rearrangeaient of terns ta eqn* CA*19) gives 

(A,f8) Slid. (A*21) are used, to rtaiise ite ladder 
stmotore* 

™ TIROmHitIT Mm 

If t Is the total time for mi» elodi oycla then p 1/f 
is the throoghoiit rat® of the stnieture* la other %«a*ds^ it is 
lisflatd as thwa Bta3i;|.a\M aimher of Mt parallel lK|>at seqaeaots ■ 
that oaa he prooessed la oae olohk oycl® hy the processing 
els«ia1»« 
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tijexi i« proceed as folloiis 

1.1 

1*2 


1.3 


Let i • o 
FiM 1^^ •» 

(I -(I 


»i » 


■ D , 
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t \*1/2 |_t 


ii " *^©i% 


-( 1 Cl 




1.4 


Obtain ^ ^ to get 



TMs $i>pp®Mix aimf& r 0 )^ 0 mm for ^Etriotis ' 

Im ia^lmmt&tiona of & given transfer fnnoti^* 

As sa ilXnstratiea^ » tliird m^er Bnttsnforth normalised 
Xo«l»sss filter is considered'* Its reference dcHsain ^tensfer 
fmmti&u is given m 

Mi^) • 

Fen* tbe direct i^eliimitiCKia Mst&eds discussed in section %* 
tJie s-«4Gffiiain transfer fimction using eqn# C3*1a) is obimined m 

E(z) w 

For quantization of tlie design psrane^rsg mmber of 
bits used are 6, 4 and 2 respectively* For ciM^parison* a 
conventional direct form implementation* using coefficients as 
multipliers* is also included# All compubsr progxnms are 
listed in Appendix C, 

As a variation* f<w? GrajHtarkel lattice filter* a third 
carder Cbebysbev nsnvdliiied lo^pass filter is consid®E'ed teving 
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c FRSSUmirY RESPONSE OF 

C CONVENTIONM, DIRECTLY IMPLEMEaWED DIGITAL FILTER 

dimension A(10),B(10), H{8, 100), FR(IOO) 
complex ZINV,NUM,DEN 
C Input Format; 

C I:flP=No, of points ,ND®3= degree of trans .fianction. 

C A(z)= nian.of trans. function; B(z )= denominator 

PI=4.*atan(l . ) . 

read {20,*} NP,NDH3 
do 1 l©T=8,2,-2 

DV=PI/(1,*NP); T=1./{1.*2**{NBT+1)) 

read {20,*),(A(T),I=1,NDEG+1); read (20,*), {B( I ) ,1=1 ,NDEG+1 } 

M’sBd) 

do 10 Isl,NDEG+l 

A( I )= ( if ix{ ( A ( I )/TMP+T ) *2 **NBT ) *1 . )/ { 1 . *2 **NBT ) 

B( I )= ( if ix( (B (I )/™p+T ) *2**NBT ) *1 . )/ ( 1 . *2 **NBT ) 

10 continue 

do 20 J=1 ,HP 

W=(J-1)*D¥ ; ZINV=cos(V)*{1.0,0.)-sin{W)*(0.,l.) 

Nl®#cfflplx{A(l ) ,0 . ) ;DEN=cmplx(B(l ) ,0 . ) 
do 30 I=2,NDEG+1 

NUM=NUM+A{I )*ZINV**{I-1 ) ;DM=DEN+B{I )*ZINV**(I-1 ) 

30 continue 

H(NBT,J}=abs(NUM/DEN) ; FR(J)=V/PI 

20 continue 

1 continue 

do 12 1=1 ,NP 

write (21,*),(FR{I),H(8,I),H(6,I).H{4,I).H{2,I)) 

12 continue 

s top 
end 

C 

C LADDER / UNIT ELEMENT BASED WAVE DIGITAL FILTER 

complex ZE,ONE,FA,FB,FC,FD,SA,DB,SC,SD 

dimension NTYPE(10 ) ,LC(10 ) ,G1{10 } ,G2{10 ) ,G3{10 ) ,NAME{10 ) 

dimension H( 10, 100), FR( 100) 

C input format ; 

C NP= no. of points ;NBT=no . of bits used;NE=no. of 'blocks ' 

C NAI® = 'ad ‘ for adapter ; 'ue ' for unit element 

C NTYPE = +1 for series adapter ; -1 for parallel adapter 

C LC = +1 for capacitor ; -1 for inductor {3rd adapter port) 

C T1,T2,T3 ; Res ./admittances for series/parallel adaptors. 

PI=4.0*atan(l,0 );ONE=(1.0,0.0) 
read (20.*) NP.NE 
do 11 NBT=8,2,-2 
TMP=1.0/2**(NBT+1) 
do 10 1=1 ,NE 

read (20,100 ) NAME(I) 

100 format(a2) 

if (NAME(I).eq. 'AD ') go to 20 
read (20 ,* ) R 
go to 10 

20 read (20.*),NTYPE(I),LC(I),T1.T2,T3 

Yl=2 .0*T1/(T1+T2+T3 );Y2=2.0*T2/{T1+T2+T3 ) 
Gl(I)=(ifix{(Yl+™p)*2**NBT)*1.0)/2**NBT^^^^^^ ^^^^^^^ 

G2{I )=(ifix( (Y2+TMP)*2**NBT)*1 .0 )/2**NBT 
G3(I)=2.0-G1(I)-G2{I) 

10 continue 



u u 
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D\f=pi/m> 

21 do 25 J=1,NP 

V^{J-1)*D¥ ; Rl=cos{W); R2=sin(W) ;ZE=cniplx{Rl ,R2 ) 
FA=ONE; EB= (0 .0 .0 .0 ) ;FD=ONE;FC=FB 
do 40 I=NE,1,-1 

if (NAME(I).eq. 'AD') go to 30 
call UNITEL{SA.SB,SC,SD,ZE) 
go to 35 

30 call ADAPTR{NTyPE(I),LC{I),Gl(I),G2(I),G 

13(I),SA,SB,SC,SD,ZE) 

35 FA=SA*FA+SB*FC;FB=SA*FB+SB*FD 

FC=SC*FA+SD*FC ; FD=SC*FA+SD*FC 

40 continue 

H( NBT , J )=ABS(ONE/FA ) ;FR{ J )=¥/PI 

25 continue 

11 continue 

do 12 1=1, NP 

WRITE {21,*),(FR(I),H(8,I),H(6,I),H(4,I).H(2,I)) 

12 continiie 
. stop ^ 

end 


FREQUENCY RESPONSE OF LATTICE WAVE DIGITAL FILTER 



100 


20 


10 


30 

35 


dimension NTyPE(10 ),LC(10 ),G1{10 ),G2{10).G3{10),NAME{10 ),FR(100 ) 
dimension HH(8,100) 

complex ZE ,ONE ,FA ,FB ,FC ,FD ,SA ,SB ,SC ,SD ,H( 2 , 100 ) 

PI=4.0*atan(1.0 ); ONE= (1 .0 ,0 .0 ) 

read (20,*) NP 

D1»P=PI/(NP*1.0 ) 

do 1 NBT=8,2,-2 

TMP=1.0/2**(NBT+1) 

do 50 LOOP=l,2 

read (20,*) NE,NT,NM 
do 10 1=1 ,NE 

read (20,100) NAME(I) 
forniat(a2) 

if (NAME(I).eq. 'AD ') go to 20 
read ( 20 , * ) R 
go to 10 

read ' ( 20 , * ) ,NTYPE( I ) ,LC( I ) ,T1 ,T2 ,T3 
Y1=2.0*T1/(T1+T2+T3) ; Y2=2 .0*T2/ (T1+T2+T3 ) 

G1 ( I )= ( if ix( ( Y1 +TMP ) *2**NBT )*1 . 0 )/2**NBT 
G2 ( I )= ( if ix ( ( Y2 +TMP ) *2 **NBT ) *1 . 0 )/ 2**NBT 
G3(I)=2.0-G1(I)-G2(I) 

continue 
do 25 J=1 ,NP 

W=1.0*(J)*DW 

Rl=cos(¥); R2=sin(W) ; ZE=cmplx{Rl ,R2 ) 

FA=ONE; FB=(0.0,0.0) ; FD=ONE; FC=EB 

do 40 I=NE,1,-1 

if (NAME(I) .eq. 'AD ') go to 30 
call UNITE(SA,SB,SC,SD,ZE) 
go to 35 

call ADAPTR(NTYPE(I),LC(I),G1(I) 

1 ;G2(I),G3(I),SA.SB,SC,SD,ZE) 

FA=SA*FAfSB*FG ;FB=SA*FB+SB*FD 
FC=SC*FA+SD*rc ;FD=SC*FB+SD*FD 
continiie 

if (NT.eq.l) go to 22 


40 



H {LOOP , J )= ( PC+FD*1JM )/ (FA+FB *NM ) 
go to 25 

22 H{L00P,J)=(FC+H)*NM/ZE)/(FA+FB*NM/ZE) 

25 continue 

50 continue 

do 27 J=1,NP 

E1?(J}=J*1,0/NP ; HH(NBT,J)=abs(H(2,J)-H(l,J)) 

27 continue 

1 continue 

do 12 1=1 ,NP 

write (21,*),(FR{I),(HH(NBT,I),NBT=8,2,-2)) 

12 continue , 

stop 
end 

C 

SUBROUTINE ADAPTRCNT.K, XI ,X2 ,X3 ,A,B ,C ,D ,Z ) 
complex DN ,C1 ,02,03 .04 , A ,B ,0 ,D , Z 
DN=1 .0^NT*K*(1 .0-X3 }/2 

01=NT*{1.0-X1)+K*(1.0-X3)/Z ; C2=-X1*{NT*1 .0-K/Z ) 

O3=-X2*(NT*1.0-K/Z) ; 04=NT*(1 ,0-X2 ) +K*{l .0-X3 )/Z 
A=DN/03 ;B=-C4/03 ,-O=01/C3 ,-D= (02 *03-01*04 )/ (DN*03 ) 
return 
end 

SUBROUTINE UNITEL (A, B,0,D,Z) 
complex A,B,0,D,2 

A=sqrt(Z);D=ONE/A ; B=(0 .0 .0 .0 ) ;0=B 

return 

end 


{ GRAY AND MARKED ALGORITHM FOR LATTIOE DIGITAL FILTER 

PROGRAM GRAYMARKEL ( INPUT ,OUTPUT ) ; 
var 

M,N,I ,J;integer ; 

TMP:real; 

K,V:array[0 . .10] of real; 

A.P,ZB;array[0..10,0..10] of real; 
begin 

readLN M; 
for N;=0 to M do 
read P[M,N]; 
readLN; 

for N;=0 to M do 
readA[M,N]; 
readLN ; 

TMP;=A[M,0]; 
for N:=0 to M do 
begin 

A[M,N]:=A[M,N]/TMP; P[M,N] :=P[M,N]/TMP 

end; 

for I:=M downto 1 do 
begin 

for J;=0 to I do ZBCX,J] •=A[I ,I-J] ; 

K[I-1];=A[I,I]; 

for J:=0 to I-l do 

A[I-1 ,J] :={A[I ,J]-K[I-1]*ZB[I .J] )/(1.0-Kf:i-13*K|I-13 ) 
V[I]:=PCI,I3; 

for J;=0 to I-l do P[I-I,J];=P[I,J3-ZB[I,J3*V[I] 

• end; 




V[0];=P[0 ,0] j 

VKITEUi? { 'K- parameters are'); 

for N;=0 to M~1 do WRITE { ' k( ’ )= ',K[N] :10 :6 ) ; 

WRITELN; 

WRITEUM( Tap parameters are'); 
for N:=0 to M do 

\®ITE( ' v( '.N;!, ')= ',V[N]:10;6) 

end. 

C — — : 

C FRl^UEIKinf RESPONSE OF GRAY AND l^aRKEL LATTICE DIGITAL FILTER 

dimension A{10 ),ZBM{10 ,10 ) ,H(8 ,100 ),FR(100 ) ,RM{10 ) ,V{10 ) 

complex ZINV,NW,Diajl,RM 

PI-4. *atan(l .); read (20,*) NP,NDEX3 

do 1 NBTs8,2,-2 

DWI/(1.*NP) 

T=1./{1.*2**(NBT+1)) 

read (20,*),(A(I),I=1,NDEG+1) ; read (20 ,*) , (V( J) ,J=1 ,NDEG+1 ) 
read (20,*). ( (ZBM(J,I ) ,1=1 ,NDEG+1 ),J=1,NDEG+1 ) 

TMP=A(1) 

do 10 I=1,NDEG+1 

do 15 J=1,NDEG+1 

ZBM( J , I ) = ( i f ix ( ( ZBM( J , I )/']MP +T ) *2 * *NBT ) *1 . ) / ( 1 . *2 * *NBT ) 

15 continue 

A(I)=(ifix((A(I)/TMP+T)*2**NBT)*l. )/(l.*2**NBT) 

10 continue 

do 20 JJ=1,NP 

W=(JJ-1)*DW ; ZINV=cos(W)*(1.0,0.)-sin(W)*(0,,l.) 
DEN=cmplx(A(l),0. ) 
do 16 J=1,NDEG+1 
RM( J )=CMPLX( ZBM( J , 1 ) , 0 . 0 ) 

16 continue 

do 30 I=2,NDEG+1 
do 17 J=1,NDEG+1 

RM(J)=RM(J)+ZBM(J,I)*ZINV**(I-1) 

17 continue 

DEN=DENfA(I )*ZINV**(I-1 ) 

30 continue 

NUM=(0.0,0.0 ) 
do 18 J=1,NDEG+1 
NUM=NUM+-RM( J )*V( J ) 

18 continue 

H(NBT,JJ)=abs(NUM/DEN) ; FR(JJ)=W/PI 

20 continue 

1 continue 

do 12 1=1, NP 

write (21,*),(FR(I),H(8,I),H(6,I),H(4,I),H(2,I)) 

12 continue 

stop 
end 

C — — — — 

C FREQUaCY RESPONSE OF LBR 2-PAIR BASED DIGITAL FILTER 

complex ZE , ONE .FA,FB ,FC ,FD ,SA,DB ,SG,SD 

dimension NAME(10 ) ,SIG(10 ) ,ALP(10 ) ,K(10 ) ,H(10 .100 ) ,FR(100 ) 
PI=4.0*atan(1.0) ; ONE=(l .0 ,0 ,0 ) 
read (20,*) NP.NE.R 
do 11 NBT=8.2,-2 

TMP=1 .0/2**(NBT+l ) ;TR= ( IFIX( (R+TMP )*2**NBT )*1 .0 )/ (2**NBT ) 
do 10 I=1.NE 



read (20,100) I) 

100 f oEwit ( a3 ) 

if (NAME( I ) .eq. 'TWO • ) go to 20 
read (20,*) K(I),SIG(I) 
go to 5 

20 read (20 .* ) ,K(I ) ,SIG(I ) ,MP(I ) 

ALP { I )= ( if ix ( ( MJP( I ) +TMP )*2 **1)3BT ) *1 . 0 )/2 **NBT 

5 SIG( I )= ( if ix{ ( SIG( I ) +TMP ) *2**NBT ) *1 . 0 )/2 **NBT 

10 continue 
D¥=PI/NP 

do 25 jn=l.NP 

^(J-1)*D¥ 

R1=co 5(W); R2=sin(W) 

ZS'caglxCRl ,R2 ) 

FA=0NE:; FB=(0.0,0.0) 

FD=ONE,v FC=(0.0,0.0) 
do 40 I=NE,1,-1 

if (NAWE(I).eg. 'TWO') go to 30 

call LBRONE(K(I),SIG(I),SA,SB,SC,SD.ZE) 

go to 35 

30 call LBRTWO(K{I),SIG(I).ALP(I),SA,SB,SC.SD,ZE) 

35 FA=SA*FA+SB*FC; FB=SA*FB+SB*FD; FC=SC*FA4-SD*FC 

FD=SC*FB+SD*FD 

40 continue 

H(NBT,J)=ABS((FC+FD*TR)/(FA+FB*TR)); FR(J)=W/PI 

25 continue 

11 continue 

do 12 1=1, NP 

WRITE {21,*),(FR(I),H(8,I),H(6,I),H(4,I),H(2.I)) 

12 continue 
stop 
end 

SUBROUTINE LBRONE{L,SG, A, B,C,D,Z ) 
complex A,B,C,D,Z 
A=(1.0+SG/Z)/1.0 
B=L*(1.0~SG)/(Z*1.0) 

C=L*(1.0-SG)/1.0 

D=(SG+1.0/Z)/1.0 

return 

end 

SUBROUTINE LBRTVD(L,SG,BET,A,B,C,D,Z) 
complex A,B,C,D,Z 

A=(l .0+BET*(l .0+SG)/Z+SG/(Z*Z ) )/l .0 
B=L*(1.0-SG)*(BET*1.0+1.0/Z)/(Z*1.0) 

C=L*(1 .0-SG)*(l .0+BET/Z )/l .0 

D=(SG*1.0+BET*{1.0+SG)/Zfl.0/{Z*Z))/1.0 

return 

end 


{ DESIGN OF ORTHOGONAL DIGITAL FILTER } 

PROGRAM VLSI (INPUT, OUTPUT); 
const 
U=2; 

type 

Y=array[0 . .2 ,0 . .5] of real; 

var 
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Kl,K2.n®>l,TMP2;real; 

T,LM©DA.TI,,TTL;Y; 

procedure MATPRO{M,P,N:integer ;var A,B,C:Y); 
var 

I ,J,K ; integer ; 
begin 

for I;=0 to M do 

for J:=0 to N do 
begin 

e[I,J3;=0.0; 
for K:=0 to P do 

C[I,J]:=C[I,J3+A[I,K]*B[K,J] 

end 

end; 

begin 

readM NDBG; 
for I :=0 to 2 do 
begin 

for J:= 0 to NDEG do 

read LAMBDA[I,J]; 

readO;! 

end ; 

WRITELN; 

MJITELif; 

’WlITEIiJ; 

WRITE { 'Coefficient values for Pipelined Orthogonal '); 

WRITELN ( 'Digital Filter : Degree = ‘,NDEG;1); 

for I;=0 to NDEG do 

begin 

K1 : =LAMBDA[ 1 , 0 3/LAMBDA[ 0 , 0 ] ; 
if Kl>1.0 then Kl:=1.0; 

If KK-l.O then Kl;=-1.0; 

TMPl : = ABS { SQRT ( 1 . 0-Kl *K1 ) ) ; 

T[0 ,0 ] :=1 .0/TMPl ;T[0 ,1 ] :=-Kl/TMPl ;T[0 ,2 ] :=0 .0 ; 
T[1,0];=-K1/TMP1;T[1,13;=1.0/TMP1;T[1,2]:=0.0; 
T[2,03:=0.0;T[2,1]:=0.0;T[2,2]:=1.0; 

MATPRO ( U , U , NDEG ,T , LAMBDA , TL ) ; 

K2;=TL[2,0i/TL[0,0]; 
if K2>1.0 then K2;=1.0; 
if K2<-1.0 then K2;=-1.0; 

TMP2 : = ABS ( SQRT { 1 . 0-K2 *K2 ) ) ; 

T[0 .0 ] ;=1 .0/TMP2 ;T[0 .1] :=0 .0 ;T[0 ,2] :=-K2/TMP2 ; 
T[1,0]:=0.0;T[1,1];=1.0;T[1,2 3;=0.0; 

T[2,0 3: = -K2/TMP2;T[2,13;=0.0;T[2,2 3:=1.0/TMP2,- 
MATPRO ( U , U ,NDEG ,T ,TL ,TTL ) ; 

WRITELN; 

WRITELN ( 'K( ',1:1, ',1 ) = '.Kl;7:4,' K{ ',1:1, ',2 ) - ',K2:7:4J; 

for J;=0 to NDES3 do 

LAMBDA[0 ,J3 ;=TTL[0 ,J3 ; 
for J;=0 to NDEG-1 do 
begin 

LAMBDA[ 1 , J 3 ; =TTL[ 1 , J+1 3 ; 

LAMBDA[ 2 . J 3 :=TTL[ 2 , J+1 3 

end; 

LA1®DA[1.NDK33:=0.0; 

LAPBDA[2,NDEG3:=0.0 

end 


end. 
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